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Abstract 


: 


The  problem  of  laser  cooling  of  a  single  ion  or  atom  in  a  harmonic  trap  was 
considered.  A  simple  theory  of  sideband  cooling  has  been  developed.  In  the  limit 
that  the  particle's  secular  motion  can  be  treated  semiclassically.  the  theory  allows 
the  calculation  of  a  logarithmic  cooling  rate  either  numerically  or.  for  two  special 
cases,  analytically.  The  theory  could  be  used  to  optimize  the  parameters  of  a 
cooling  experiment. 

The  spectroscopy  of  a  single  atom  in  the  ladder  configuration  has  been 
treated  theoretically.  A  dressed  atom  approach  was  used  to  provide  qualitative 
information  about  the  system.  The  optical  Bloch  equations  for  the  four  level  sys¬ 
tem  in  the  rotating  wave  approximation  were  developed  and  solved  for  steady 
state.  The  Bloch  equations  were  also  solved  in  the  adiabatic  approximation  and 
upward  and  downward  transition  rates  were  extracted  from  this  treatment. 


I.  Introduction 


A.  Motivation 

Breakthroughs  in  the  trapping  and  cooling  of  ions  have  made  it  possible  to 
perform  spectroscopic  experiments  on  a  single  ion  trapped  indefinitely  in  a  volume 
with  fim  dimensions  (31).  Recent  developments  in  the  trapping  and  cooling  of 
clouds  of  neutral  atoms  (26)  indicate  that  confining  a  single  neutral  atom  is  feasi¬ 
ble.  Trapped,  cooled  atoms  will  possibly  provide  the  ultimate  laser  frequency 
standard.  Alone  in  the  trap,  the  spectral  lines  of  the  atom  will  not  be  collisionally 
broadened  and  if  restricted  to  a  region  smaller  than  X/4,  the  Doppler  effect  will  be 
suppressed  (12).  It  is  also  possible,  using  the  "atom  amplifier"  technique  first  pro¬ 
posed  by  Dehmelt  (13),  to  use  extremely  weak  transitions,  accessing  states  with 
extremely  long  lifetimes  and  consequently  narrow  natural  linewidths. 

The  theory  of  laser  cooling  in  a  harmonic  trap  has  been  developed  by  several 
researchers  using  different  methods  (10),  (11),  (18),  (19),  (20),  (2 1 ).  (29).  (36).  (38). 
(39).  (-40).  Most  of  these  theories  are  cumbersome  and  involved.  A  simple  theory 
that  displays  the  physical  mechanisms  of  laser  cooling  and  is  numerically  correct 
in  regions  of  interest  is  desirable. 

It  is  possible  to  observe  effects  in  the  spectroscopy  of  a  single  atom  that  are 
masked  when  the  spectroscopic  sample  contains  many  atoms.  One  of  the  most 
interesting  examples  of  this  which  has  been  observed  is  the  phenomenon  of  "elec¬ 
tron  shelving"  (3),  (26).  (33).  In  an  electron  shelving  experiment  (also  known  as 


quantum  jump.  A  simple  example  of  such  an  experiment  uses  the  ’V" 
configuation.  shown  in  Figure  1.  Two  lasers  are  focused  on  the  atom,  one  exciting 
a  strong  transition  (Einstein  coefficient  .-i^lOss-1),  and  one  exciting  a  weak  transi¬ 
tion  (Einstein  coefficient  A  ^=ls-1).  The  transitions  share  a  common  lower  level. 
Ordinarily  the  strong  transition  scatters  about  10s  photons  per  second,  which  can 
be  observed  with  a  photodetector.  If.  however,  the  atom  makes  the  transition  to 
the  upper  state  of  the  weak  transition  (state  2),  then  the  atom  would  stop  scatter¬ 
ing  photons  on  the  strong  transition  until  the  atom  decayed  back  into  the  lower 
state  (state  0).  Thus  single  quantum  jumps  on  a  very  weak  transition  can  be 
observed.  Other  configurations  of  the  atoms  energy  level  structure  have  been  con¬ 
sidered  theoretically.  Some  examples  are  the  "A"  configuration  which  is  an 
inverted  "V"  (22),  (32)  and  various  three-level  cascades  (32). 

One  configuration  which  has  not  previously  been  considered,  but  which  pro¬ 
vides  some  interesting  details  is  a  four-level  cascade  in  which  two  strong  transi¬ 
tions  are  coupled  by  a  weak  transition.  This  ladder  configuration  is  illustrated  in 
Figure  2.  One  of  the  debates  about  the  experiments  involving  three-level  electron 
shelving  schemes  such  as  the  ’V"  configuration  has  to  do  with  quantum  measure¬ 
ment  theory.  Currently  accepted  interpretations  of  quantum  theory  hold  that 
quantum  systems  do  not  exist  in  any  particular  quantum  state  until  a  measure¬ 
ment  projects  them  into  some  state.  The  question  is  whether  a  period  of  dark¬ 
ness,  the  absence  of  a  fluorescence  signal,  constitutes  a  measurement  which  pro¬ 
jects  the  atom  into  the  metastable  state,  state  2  of  Figure  1.  For  the  ladder 
configuration  there  are  no  periods  of  darkness,  instead  the  fluorescence  alternates 
between  wavelength  -y  when  the  weak  transition  is  not  excited,  and  wavelength 
-*-3  when  the  weak  transition  is  excited.  A  change  in  the  wavelength  of  the  fluores¬ 
cence.  which  is  a  positive  signal,  as  opposed  to  the  absence  of  a  signal  for  three- 
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indicates  a  quantum  jump.  Thus  an  experiment  '.using  the  laii-r 
configuration  could  resolve  whether  what  is  being  observed  is  indeed  a  single 
quantum  event. 


B.  Problem  Statement  and  Approach 

The  purpose  of  this  research  is  twofold.  First,  in  Chapter  III.  a  simple  theory 
of  sideband  coding,  in  the  limit  that  the  amplitude  of  the  ion’s  motion  is  less 
than  the  wavelength  of  the  cooling  light  (the  Lamb-Dicke  limit),  is  developed. 
This  is  accomplished  by  first  defining  the  operators  which  describe  a  two  level 
atom  in  a  harmonic  trap  and  writing  the  Heisenberg  equations  of  motion  for  these 
operators.  The  underlying  assumption  that  makes  it  possible  to  simplify  the 
theory  is  then  made.  This  assumption  is  that  the  atom  is  "well  localized",  i.e..  the 
uncertainty  in  the  expectation  value  of  the  operator  corresponding  to  the 
particle's  position  vector.  T-.  is  small  compared  to  the  wavelength  of  the  cooling 
laser,  so  that  r*  can  be  replaced  by  its  expectation  value,  Cr*>.  when  taking 
expectation  values  of  operators  that  are  functions  of  T.  The  equations  of  motion 
for  the  expectation  values  of  all  the  pertinent  operators  then  follow  directly  from 
the  Heisenberg  equations  of  motion  for  the  operators.  It  is  then  possible  to  derive 
three  simple,  first  order  Bloch  equations  for  the  twolevel  atom  in  a  harmonic 
trap.  For  one  dimensional  cooling,  a  dimensionless  displacement  described  by  a 
second  order  differential  equation  is  defined.  The  problem  is  thus  reduced  to  a  set 
of  four  coupled,  nonlinear  differential  equations,  one  of  which  is  second  order. 


For  two  special  cases  these  equations  are  solved  analytically, 
the  weak  field  ease,  in  which  the  cooling  radiation  is  taken  to  be  s 
intensity  that  the  probability  that  the  atom  is  in  the  ground  is  v 
This  eliminates  one  of  the  four  coupled  differential  equations 
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Hgure  2.  Ladder  Configuration. 


complexity  of  the  remaining  three.  An  analytical  solution  can  then  be  derived. 

The  second  case  solved  analytically  makes  use  of  an  adiabatic  approximation. 
Such  an  approximation  is  useful  when  trying  to  solve  a  coupled  set  of  differential 
equations  that  can  be  divided  into  two  groups.  The  first  group  consists  of  equa¬ 
tions  that  have  "large"  self  relaxation  terms.  The  equations  in  the  second  group 
represent  variables  that  change  more  slowly  than  the  var'abies  described  by  the 
equations  in  group  one.  If  one  is  interested  only  in  processes  having  time  constants 
much  smaller  than  the  coefficients  of  the  large  seif  relaxation  terms,  then  it  is  an 
excellent  approximation  to  assume  that  the  first  group  of  variables  adiabatically 
follow  the  mere  slowly  varying  second  group.  One  can  then  use  the  steady  state 
values  of  the  first  group  of  variables  in  the  equations  for  the  mere  slowly  varying 
second  group.  This  once  again  reduces  the  number  and  complexity  of  equations. 
For  both  special  cases  considered,  it  is  shown  that  the  amplitude  of  the  dimen¬ 
sionless  displacement  decays  exponentially  and  an  analytical  expression  for  the 
decay  constant  is  presented.  Decay  of  this  amplitude  corresponds  to  decreasing 
the  atoms  translational  energy,  which  is  what  is  meant  by  cooling. 


The  four  original  equations  are  solved  numerically  using  a  fourth  order 
Runge-Kutta  method,  and  the  validity  of  each  of  the  two  analytical  models  is 
tested  by  comparing  it  to  the  numerical  results.  Exponential  decay  is  observed  for 
all  the  numerical  cases  considered,  and  a  decay  constant  is  calculated.  Good 
agreement  between  the  numerical  and  analytical  results  is  seen  in  regions  where 
the  approximations  used  in  deriving  the  analytical  results  are  valid.  Plots  of  the 
decay  constant  as  a  function  of  the  various  parameters  of  the  four  coupled  equa¬ 
tions  can  be  easily  generated,  and  some  representative  examples  of  such  plots  are 
presented. 


In  the  second  part  of  the  research,  contained  in  Chapter  IV.  the  fluorescence 
signal  of  a  single  atom  for  the  ladder  configuration  is  analyzed  theoretically.  It  is 
expected  that  the  fluorescence  will  alternate  between  two  wavelengths.  and  ~3. 
A  spi  itting  of  the  weak  transition  similar  to  the  Autler-Townes  doublet  seen  in 
the  "V"  configuration  (23)  is  expected.  A  dressed  atom  approach  is  used  to  obtain 
qualitative  information  about  this  splitting. 


The  statistics  of  fluorescence  are  determined  from  deriving  optical  Bloch 
equations  for  the  four-level  system.  The  equations  of  motion  for  the  density- 
matrix  of  a  coherently  excited,  multi-level  system,  in  the  electric  dipole  approxi¬ 
mation  are  taken  from  the  literature  (27).  The  on  diagonal  elements  of  the  den¬ 
sity  matrix.  pnn.  are  just  the  probabilities  to  be  in  the  state  n.  Only  three  cf  four 
probabilities  are  independent  since  the  total  probability  must  equal  unity.  The  off 
diagonal  elements  of  the  density  matrix,  the  coherences,  oscillate  at  the  transition 
frequencies  between  the  two  states  that  they  correspond  to.  in  this  case  optical 
frequencies.  Since  oscillations  at  such  high  frequencies  have  little  effect  on  the 
probabilities,  slowly  varying  coherences  are  defined  so  that  the  rotating  wave 
approximation  can  be  made.  In  the  rotating  wave  approximation  functions  oscil¬ 
lating  at  optical  frequencies  are  replaced  by  their  time  averaged  values. 

The  resulting  equations  of  motion  are  solved  first  in  steady  state.  The  steady 
state  solution  is  compared  with  the  qualitive  results  obtained  from  the  dressed 
atom  approach.  Finally,  the  Bloch  equations  for  the  four-level  atom  nr-  solved  in 
an  adiabatic  approximation.  The  adiabatic  approximation  is  appropriate  here 
because  the  Einstein  A  coefficients  appear  in  self  relaxation  terms  f  r  Boh  tin 
probabilities  and  the  slowly  varying  coherences.  We  have  specifically  oh  >eti  tran¬ 
sitions  so  *hat  two  of  the  A  coefficients  involved  are  very  much  larger  than  tin. 
third.  It  is  shown  that  eventuallv  all  of  the  variable-  .\un  !•<• 
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adiabaticaily  except  one.  The  remaining  variable  corresponds  to  an  inversion  on 
an  effective  two  level  system  driven  by  a  rate  process.  The  theory  of  two  level 
systems  driven  by  rate  processes  is  well  understood  (9),  (23),  and  the  statistics  of 
fluorescence  for  the  ladder  configuration  follow  directly  from  previous  work  (9). 
(23). 


13.  Background 


A.  Laser  Cooling 

1.  Experiments 

In  1980  Neuhauser.  et  al  reported  trapping  a  single  barium  ion  in  a  Paul 
radio  frequency  quadrupole  trap  using  laser  cooling,  and  holding  it  in  a  volume 
with  dimensions  of  (31).  The  Paul  trap  consists  of  three  electrodes:  two  cap 
electrodes,  shaped  like  hyperbolas,  with  a  ring  electrode  between  them.  An  oscil¬ 
lating  electric  field  is  applied  to  the  electrodes.  The  trap  is  illustrated  in  Figure  3. 
The  potential  created  in  the  trap  space  is 

V{i?,t)=v{J?)cosP.t=VQ(X'2+Y2-2Z2)cosnrf  t  (2.1) 

where  V  is  the  electric  potential.  V'0  is  a  constant,  and  Qr/  is  the  frequency  of 


oscillation  of  the  electric  field  (11).  The  origin  of  coordinates  is  at  the  center  of 
the  ring  electrode  with  the  z-axis  running  between  the  two  cap  electrodes,  perpen¬ 
dicular  to  the  plane  of  the  ring  electrode. 

A  Paul  trap,  which  has  an  oscillating  electric  field,  was  used  instead  of  a  trap 
with  a  static  electric  field,  such  as  a  Penning  trap,  because  static  traps  must  use 
both  an  electric  and  magnetic  field.  The  magnetic  field  causes  a  Zeeman  ►-•fiPct 
which  is  undesirable  for  high  resolution  spectroscopy  (11).  It  has  been  shown  (10» 
that  the  time  varying  potential  of  the  Paul  trap  is  equivalent  to  a  time  in  depen- 
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dent  effective  potential  given  by 


Vtff(R)=—±r(X-+^Zi) 


which  is  a  harmonic  potential.  The  angular  frequencies  of  oscillation  in  the  trap 
are  given  bv 


Uz  =U,J  =y2  V.  =- 


In  the  trap  used  by  Xeuhauser.  et  al  the  spacing  between  the  two  cap  elec¬ 
trodes  was  0.5  mm.  and  the  inner  diameter  of  the  ring  electrode  was  0.75  mm. 
The  frequency  of  oscillation  of  the  electric  field,  .  was  2rr  x  ISMHz  and  the 
angular  frequency  of  oscillation  of  the  ion  in  the  well,  in  the  X  and  Y  directions. 
u;,  was  2t  x  2.4\IHz.  A  barium  atomic  beam  was  injected  into  the  trap  from  an 
oven  and  ionized  by  an  electron  beam.  The  ion  was  thermalized  by  allowing 
helium  buffer  gas  into  the  trap  at  a  background  pressure  of  dCr4  Torr  (30).  A 
laser  tuned  to  the  Ba+  line  at  493.4  nm  was  shone  on  the  trap,  and  when  a 
barium  ion  was  present  in  the  trap  it  could  be  observed  visually  by  fluorescence 


The  493.4  nm  line  for  Ba~  is  a  transition  between  the  ground  6"5.,  state  and 
an  excited  6:/\  state.  For  the  purposes  of  cooling  the  ion.  the  laser  was  tuned  to 
a  frequency  below  the  resonant  frequency  of  the  transition  by  an  amount  v.  One- 
third  of  the  time,  the  ion  spontaneously  decayed  to  a  5 ~D3  .,  met ast able  excited 
state  instead  cf  b-mk  to  the  ground  state.  A  second  laser,  tuned  t he  XI)  > 
transition  at  649.9  nm  kept  the  ion  from  sitting  in  the  metastabl--  svs'c.  A n 
energy  level  diagram  is  shown  in  Figure  4.  Both  laser  beams  entered  ;!;••  tra.p 
between  the  ring  e!*-etrode  and  one  of  the  cap  electrodes,  at  an  angle  of  4-V  to  lh- 
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2.  Laser  Cooling  Theories 


The  idea  of  using  laser  light  tuned  below  resonance  to  cool  atoms  was  nrs: 


proposed  by  Hansch  and  Schawlow  (17)  for  free  atoms,  and  by  \ Vineland  ant 


Dehmelt  (3S)  for  trapped  ions.  There  is  a  fundamental  difference  between  tht 


laser  cooling  of  free  atoms  and  trapped  atoms,  at  least  in  the  Lamb-Dicke  limit. 


Laser  cooling  of  a  free  atom  takes  place  because  of  the  atom's  Doppler  shift  (171. 


For  a  running  wave  laser  field,  when  the  laser  is  tuned  to  the  lower  half  of  the 


atom's  Doppler  spectrum,  the  atom  is  shifted  into  resonance  when  it  moves 


toward  the  laser  and  out  of  resonance  when  it  moves  awav  from  the  laser.  If  the 


itom  radiates  spontaneously  in  three  dimensions,  then  the  net  effect  of  absorbing 


photons  as  it  moves  toward  the  laser  is  to  slow  the  atom,  while  the  net  effect  of 


absorbing  photons  as  it  moves  away  from  the  laser  is  to  speed  the  atom.  Because 


of  the  Doppler  shift,  the  atom  absorbs  more  photons  when  it  is  moving  toward 


the  laser  than  awav  from  it.  so  the  atom  is  slowed.  Laser  cooling  of  free  atoms 


and  atoms  in  radiation  traps  has  been  treated  extensively  in  the  literature  (2).  '5'. 


(6).  (7).  (S).  (16).  (25),  and  will  not  be  discussed  further  here,  except  to  note  that 


in  reference  (3)  Cook  uses  Ehrenfest’s  theorem  and  optica!  Bloch  equations  to  .leal 


vith  atomic  motion  in  resonant  radiation.  The  method  used  to  come  up  with  a 


simple  theory  of  side-band  cooling  used  in  this  work  can  be  thought  of  as  an 


extension  of  Cock's  method  to  bound  atoms. 


A  Doppler  shift  picture  is  a  good  physical  description  cf  the  cmTing  of 


trappe  1  a*  un  as  long  as  the  volume  to  which  the  atom  is  con!: nod  is  larg- 


pai'e- i  t  ■  th“  waveh-ng* h.  X.  of  the  cooling  laser.  However,  wlvn  an 


is  coni: i  r  >  a  vo’.um**  whoso  dimensions  are  smaller  than  X  1  D  shift  h 


uppro-s>:-d  (12!.  so  las^r  cooling  of  trapped  ions  in  this  limit,  the  Lamb-Dick- 


e  ex;  lain-!  by  the  D  ppl-r  shift.  Laser  r  .-cling  in 
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Lamb-Dicke  limit  takes  place  via  side-band  ceding.  The  electronic  levels  of  the  icn 
are  split  by  the  energy  levels  in  the  trap.  For  a  harmonic  trap  these  levels  are 
evenly  spaced  tu  apart,  where  w  is  the  angular  frequency  cf  oscillation  of  the 
atom  in  the  trap.  Figure  5  is  an  energy  level  diagram  for  a  two  state  atom  split 
by  the  trap  levels.  The  spectrum  of  the  ion  then  looks  like  Figure  6.  with  the 
central  peak  at  ~0  and  sidebands  at  ~0 ±m  u  where  m  is  an  integer.  When  the 
laser  is  tuned  to  a  lower  sideband  (w0— mv\  the  ion  is  excited  from  the  | j.n> 
state  (ground  electronic  state,  nth  harmonic  oscillator  level)  to  the  |  e.n—m> 
state  (excited  electronic  state,  n—m  harmonic  oscillator  level)  from  which  it 
decays  to  \g.n—  m>.  Thus,  the  net  effect  is  to  decrease  the  expectation  value  of 
the  harmonic  oscillator  number  operator.  <fn>.  which  corresponds  to  cooling  the 
ion's  secular  motion. 

The  theory  of  cooling  bound  atoms,  and  specifically  harmonically  bound 
atoms  has  been  treated  by  several  groups.  Stenholm,  Javanainen.  and  Lindberg 
(19).  (20),  (21)  have  worked  in  this  area  both  together  and  separately  for  many 
years.  Reference  (21)  used  a  density  matrix  approach  to  deal  with  the  final  stages 


of  laser  cooling  of  harmonically  trapped  ions.  This  approach  involved  coupled 
equations  for  all  the  harmonic  oscillator  states  which  even  in  the  Lamb-Dicke 
limit  can  mean  hundreds  of  equations.  These  equations  can  be  solved  numerically 
for  particular  cases  cf  interest.  No  attempt  was  made  in  this  approach  to  deter¬ 
mine  a  cooling  rate.  The  parameter  of  interest  to  these  authors  was  the  minimum 
temperature  achievable  at  the  limits  of  cooling,  which  was  plotted  as  a  function  of 
the  detuning  for  one  set  cf  parameters  (21).  The  minimum  temperature  ’-vus 
achieved  for  the  detuning  approximately  equal  to  minus  the  angular  frequency  of 
oscillation  in  the  trap.  Stenholm  has  compiled  a  review  of  laser  cooling  theory  (29  i 
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Wiaeland  and  Itano  (IS!,  i'29).  (39)  and.  more  recently.  \ Vineland.  Itano. 
Bergquist.  and  Huiet  (40)  have  used  both  energy  rate  equations  and  quantum- 
mechanical  perturbation  theory  to  calculate  minimum  energies  ci’  ceding  and 
dE d t.  where  E,  is  the  energy  of  the  secular  motion  along  one  of  the  harmed: 
oscillator  axes.  In  references  (IS)  and  (39).  the  limiting  assumption  is  made  that 
the  cooling  transition  is  not  saturated.  These  authors,  like  Stenholm.  et  a!,  are 
chiedy  concerned  with  the  minimum  temperature  achievable  by  laser  coding,  arm: 
they  do  not  calculate  or  plot  rates  of  cooling  as  functions  of  the  various  parame¬ 
ters  of  the  problem.  However,  where  possible,  the  functional  form  of  their  expres¬ 
sion  for  dE ;dt  will  be  compared  to  the  cooling  rates  derived  in  Chapter  III. 
References  (39)  and  (40)  deal  with  the  limits  of  cooling  in  both  radio-frequency 
quadrupole  traps  and  Penning  traps.  The  authors  conclude  that  cooling  to  the 
zero  point  energy  should  be  possible,  that  is.  the  expectation  value  of  the  har¬ 
monic  oscillator  number  operator,  n=<h>,  obeys  the  relationship.  n<<  1.  They 
show  quantitatively  that  it  should  be  possible  to  detect  when  the' ion  has  reached 
the  ground  state  because  the  intensity  from  the  lower  sideband  should  vanish.  A 
qualitative  argument  that  this  would  be  the  case  was  made  by  Wells  in  reference 
(37). 


B.  Single  Atom  Spectroscopy 


1.  Theory  of  Electron  Shelving 

The  idea  cf  electron  shelving  was  first  suggested  by  DHondt  in  ] 973  as  a 
scheme  for  making  more  accurate  .atomic  clocks  (13).  It  was  not  until  ]9>3  t ha‘ 
Cook  and  Kimble  (9)  first  theoretically  treated  the  statistics  of  t he  proposed 
Thev  considered  the  ’V"  configuration  shown  in  Figure  1.  The  pm:, 


experiment. 


paper  used  rate  equations  and  looxed  at  the  ease  cl'  incoherent  excitation.  This 
paper  predicted  the  fluorescence  from  the  strong  transition  -.vould  flash  on  and  or: 
producing  a  random  telegraph  signal.  The  statistics  of  the  fluorescence  such  a; 
the  distribution  of  on  and  off  times  and  the  twotime  intensity  correlation  func¬ 
tion  were  calculated.  Sehenzle.  DeYoe.  and  Brewer  (34)  used  a  quantum  statistical 
approach  and  calculated  a  second  order  correlation  function  of  the  fluorescence 
intensity  to  show  that  a  single  atom  undergoes  "large  flucuaticns  between  a  state 
of  full  emission  in  the  strong  transition  to  a  state  of  no  emission." 


Coherent  excitation  has  been  treated  by  many  authors.  Kimble.  Cook,  and 
Wells  (23)  showed  that  for  the  ’V"  configuration  when  the  strong  transition  is 
highly  saturated,  the  excitation  of  the  weak  transition  can  be  treated  as  a  rate 
process,  so  that  much  of  the  theory  of  reference  (9)  applies  to  the  coherent  case. 
Cohen-Tannoudji  and  Dalibard  (4)  also  treated  the  "V"  configuration.  Using  a 
dressed  atom  approach,  they  also  showed  that  the  strong  transition  fluorescence 
would  have  bright  and  dark  periods  and  calculated  the  statistics  of  the  fluores¬ 


cence.  Sehenzle  and  Brewer  (33)  used  photon  counting  statistics  to  arrive  at  t! 


re 


same  conclusions.  Javanainen  (22)  used  the  quantum  regression  theorem  to  calcu¬ 
late  the  photon-counting  statistics  for  the  "A"  configuration.  Pegg.  Louden,  and 
Knight  (32)  treated  not  only  the  ’V"  and  "a"  configurations  but  also  two  three- 
level  cascades,  one  with  the  upper  transition  as  the  strong  transition  and  one  with 
the  upper  transition  as  the  weak  transition.  They  arrived  at  the  erroneous  con¬ 
clusion  that  there  would  be  no  long  dark  periods.  Arecchi.  et  a.l  (1)  showed  that 
for  the  mimeses  of  a  fremmnev  standard  a  nuked  excitation  toehniom-  has  soi:m 


I 


2.  Exeperiments  That  Demonstrate  Electron  Shelving 


Experimental  verification  of  electron  shelving  has  new  been  reported  by  three 
groups  of  researchers.  Using  the  same  experimental  arrangements  as  for  the  19S0 
Xeuhauser.  et  al  paper  (31),  Xagourney,  Sandberg,  and  Dehmelt  (2S),  have 
observed  bright  and  dark  periods  in  the  fluorescence  of  a  ion.  The  energy 
level  diagram  for  this  experiment  is  shown  in  Figure  -4.  The  Snfo  is  the  ground 
state.  The  ion  was  cooled  via  the  6 2P-+  to  ground  state  transition  by  a  stabilized, 
single  frequency  dye  laser.  As  was  mentioned  in  section  II.A.l  there  is  a  one  in 
three  chance  that  the  ion  will  spontaneously  decay  to  the  metastable  5 ~D3  state, 
so  a  second  stabilized,  single  frequency  dye  laser  was  used  to  keep  this  level  clear. 
The  shelf  level  was  the  5 ~Dh  2  level  reached  via  the  6 •Pi  2  level.  The  branching 
ratio  for  decay  to  the  5 ‘Dh  »  level  from  the  6 level  was  1  3.  The  light  source 
used  to  excite  the  ion  to  the  6 rP3/2  level  was  a  barium  hollow  cathode  lamp,  an 
incoherent  source.  Figure  7  is  a  plot  of  intensity  versus  time  from  reference  (2$). 
Note  that  when  the  ion  is  in  the  shelf  state  for  this  configuration,  there  is  no 
resonant  radiation  driving  the  ion,  so  that  the  only  way  for  the  ion  to  leave  the 
shelf  state  is  by  spontaneous  emission.  The  length  of  the  dark  periods  is  distri¬ 
buted  exponentially:  making  it  possible  to  directly  measure  the  spontaneous  decay 
lifetime  of  the  5 'Dhri  state.  The  spontaneous  decay  lifetime  measured  for  the 
o2Ds  state  was  greater  than  30  seconds. 

Sauter,  et  al  (33)  have  done  a  similar  experiment  with  B;U  making  use  of  a  A 
configuration  shown  in  Figure  S.  Figure  9.  from  reference  (33\  is  a  plot  of  fluores¬ 
cence  versus  time  for  a  single  trapped  ion. 


Bergquist.  et  al  (3)  have  performed  a  very  simular  experiment  us: 


cury  ion.  Their  experiment  actually  uses  a  "V"  configuration.  Figure  10  is  a;, 
energy  lev-.-!  diagram  for  this  system.  Figure  11  is  a  plot  of  intensity  ver.~u~ 
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fluorescence  photon  counts/sec 


Figure  7 


Intermittent  !'  luorpscence  from  Nngmirney,  et  nl  (2$). 
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SINGLE  ION  FLUORESCENCE 


trcm  reference  :  3 ' .  Tne  h:et;me  c!  tne  sneiving  state  :s  aoout  0.1  secor.  i.  m 


shorter  than  for  Naacurnev.  et  al.  Fhaure  11(a)  is  a  plot 


>f  the  ffucr1 


without  the  "shelving"  loser  turned  on.  The  jumps  in  the  hourescer.ce  are  iue  to 
off  resonance  transitions  caused  by  the  "strong"  transition  laser.  Figure  li  b1  is  a 
plot  of  the  intensity  for  a  single  ion  in  the  trap.  Note  that  there  are  two  levels  of 


intensity,  which  correspond  to  the  atom  on  and  the  atom  off.  Figure  11:  ;)  is  a 
plot  of  the  intensity  for  two  trapped  ions.  Note  that  there  are  three  Tseerr.abT 
levels  of  intensity,  which  correspond  to  both  ions  on;  one  ion  on.  one  ion  off:  an  .i 
both  ions  off. 


•» 
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ID.  A  Simple  Theory  of  Sideband  Cooling 


A.  Theory 

The  operators  which  describe  a  harmonically  trapped,  cooled,  twc-leve!  ic: 
can  be  divided  into  three  groups.  First  there  are  the  operators  which  describe  th 
center  of  mass  motion  of  the  ion,  7*  and  ~p.  In  the  absence  of  any  driving  force 
the  equations  of  motion  for  these  operators  are  determined  by  the  unperturbc. 
hamiltonian  for  the  translational  motion 


r'r  1  /  .  °  °  ,,  f  o  .  o  o  .  ■■ 

— (Pz'+P-J  -P: 
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where  u. .  i /.r  and  v.  are  the  frequencies  of  oscillation  of  the  ion  in  the  x.  j.  and  . 
directions  respectively.  Xext  there  are  the  operators  which  describe  the  interna 
motion  of  the  ion.  For  an  ion  with  only  two  states.  |l>  and  1 2 > .  these  opera 


tors  mav  be  taken  to  be 


<7=  1X2  . 


a-+=  [2><i  |. 


(3.2b 


5=12X21-11X1 


.'Ov’CvW'Ofv ’■•i'iv'i- -cVl- vvih  -.wC-l  -2 


3.3 


The  hamiitonian  for  the  unrerturbed  internal  motion  reads 


3, 


where  _n  is  the  frequenev  of  the  transition  between  states  1 2 >  and  |l>.  Finallv. 


there  is  the  interaction  between  the  laser  field  and  the  ion.  A 


classical:' 


prescribed  field  is  assumed.  For  a  classically  prescribed  field,  in  the  electric  dipole 
approximation  the  interaction  hamiltonian  is 

H  ).  (3.4  ) 


Here  T.T.t  ,  is  the  applied  electric  field  and  ~p.  is  the  electric  dipole  oper: 
two-level  atom  which  mav  be  written  as 


.ter  lor  a 


?=TH&++&),  (3.5) 

where  7T=<l|/T|3>  is  the  dipole  transition  moment.  For  the  purposes  of  this 
paper  7?  is  assumed  to  be  real.  For  a  running  wave  laser  field.  can  be  taken 

to  be 


T C T,  t  )=e^-  E ccs ( k  ■?—  ~t ), 


where  TV  is  the  unit  polarization  vector.  If  n  is  defined  as  7 ?■?£,  then  H  is  gi 


H  =— /i£,(d~-rd)cos(\--7*— -•<). 


The  total  Hamiltonian  of  the  svstem.  //.  is 


H  =Ht  +H,+H  . 


The  Heisenoerg  equations  of  motion  ior  the  operators  7*.  ~p.  a.  and  a3  are 

dT_  y.H  _  p* 
d  t  1 ft  m 

-^y-=— m  {v~xe,  —u'ye.j  u:zT.)—k  jJ.E(c 7+-Hj)siru  k  -T—~t ). 


da  .  -v-E  . 

— =-1  ^0a-i-—a3cos(k  ■  r -of , 
d  t  n 


—=2i-^-{a+-a)cos{T-t-^t). 
dt  n 


At  this  point  the  approximation  is  made  that  the  atom  is  well  localized,  ti 


Ar=V  \<t>-<t>-<tt>  I  (3.: 

is  small  compared  to  the  wavelength,  \=2~/k.  This  means  that  when  taki 
expectation  values  r*  can  be  replaced  by  <r*>  in  si n(V-T—~t)  and  cos(T T—~ 
The  equations  of  motion  for  the  pertinant  expectation  values  become 

4£-£.  (3.1: 

d  t  m 

—£-=—m(v;xTz  -rVy'jEy  +u:zT.  )—Vi.iE (a’  -~a)s\n(TT~^t ),  (3.1 1 


~-=-{,3+i^0)a-i-^azcos[k-  r  -~l ), 


•.’W 


- ='j j- - 1 a  —O'cos\  rc  r  i  ji o-,—  1 1. 
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'.vhere  T=<7*> .  ~p=<~p>.  o=<o>.  o3=<oz>.  Spontaneous  emission  has  net  b 
included  in  the  Hamiltonian.  Spontaneous  is  included  by  following  the  treatm 
of  reference  (‘27).  Including  spontaneous  emission  will  change  equations  (3.1 


and  (3.1  Id)  to 


do  ,  .  , , 

—  =  -(  J-i  ~'q)0  1— — — (73COSl  /v  ■  r  -af ) 

d£  7T 


d°3  0  .  /jE 

- =2: — — 

d£  o' 


(o+-<T)cos(T-T-~t')-23(o3+I) 


where  3  is  half  the  Einstein  A  coefficient. 

It  is  convenient  to  introduce  slowly  varying  Bloch  variables  u.  v.  and 
defined  as 


u=o’ e  iUJt  T-ae1'1 , 


v=—i{o,e~i'jJt  — <7e’"’( ), 


w  — o3. 


In  the  Lamb-Dicke  limit  (T 7*<<l)  we  may  expand  the  trigonometric  function: 
equations  (3.11)  and  (3.12)  to  first  order  in  the  small  quantity  V-T. 

sin (k  -T— a£  )=sin(k  -T)cos ~t  —cos [k  T)sina£  ^z(k  -T)cosaf  —sin  of ,  (3.1 


«T(TiTt'r  I'tyrr 


fft'iTWfCn.TOT.n 


The  equations  of  motion  for  the  Bloch  variables  and  p  become 


- =—3u  -f-Ar  —Clw  f  2  k  r  sin '-7  —  sin 2  A  ). 

c if  { 


1 3.15a 


— -=-Au  —it-  ~Q(k  -Tsin2~f  —  2cos"~'f ). 
at 


(3.1ob 


-=— fi  2{k-Tu  sin2  ~t  —  ccos2~f )—  ( u  —  k  •Tv  (sin  2  ~f  —  2i(  r  —1  i 


-=— m  {v:fe'z  —u2’Je'.,  ) 


(3.15.: 


-BrH  k  ■  ru  cos"_'f  —  vsin'-A  — Vi(  u  -rk  -Tr)sin2~f  , 


where  A=w— w0  is  the  detuning  and  P.=pE  It  is  the  Rabi  frequency. 


The  terms  sin'~C .  cos 2~f ,  sin2 ~t .  and  cos2 ~t  vary  much  faster  than 


n  ,  v .  w 


and  ~p  so  that  only  their  time  averaged  values  have  much  effect.  In  the  rotating 
wave  approximation  sin’-vf  and  cos ~^t  are  replaced  by  their  time  averaged  values 
l2.  Sin2~f  and  cos2.vf  average  to  0.  We  consider  cooling  of  the  x  degree  of  free¬ 
dom.  that  is  T  is  taken  in  the  x  direction  [k.=k,  k,j=kz= 0  and  i ■>x=u).  The  equa 
tions  of  motion  are  simplified  by  introducing  the  dimensionless  displacemen 
s  =kx .  The  equations  of  motion  become 


- — — 3 u  -r  At’  -f-Q il's  , 


(3.16a 


VCa  1 


w;i 


— =-A  y-dr-fiir. 
dl 


o.  i  D : 


d  w 
dl 


=  -n  o U  -v)-2  S(  T-l 


d‘5  •■>  , .  •  • 

— --r-v~s  =— A  ;  ou  — !. 

di ' 


3.1  '5 


where  K=1rk': 0  [2m).  This  set  of  four  coupled,  nonlinear  differentia! 
describe  the  cooling  process  within  the  approximations  given  above. 

The  steady  state  solution  of  equations  (3.16).  u0.  r 0.  tr0,  and  s0.  :c 
(L.amb-Dicke  limit)  is 


UQ  = 


V0~ 


v  0=- 


fX^+.d*  c 


d2+A2+H>n2 

Odd-Asp) 
d-+A2-rV'2p.2  ' 

(d--A2) 

"  d2+A2-fV2n- 


(3.17 


(3.17b 


(3.17 


5q- 


nk-3n- 


2m  id" ( d * + A" -f  - V2 IT 


(3.i; 


Other  solutions  to  equations  (3.16)  can  be  written  in  the  form 


f(t)=/o+ft(n 


(3.1  S) 


•’nerrt  /.!?  )  has  a  steady  value  of  zero.  When  solutions  of  this  form  are 
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* 


tutea  into  equation  (3.1b)  we  set 


d  u,  _ 

— — =— iu,  ~Av,  —  flsgir.  — 
dt  ’ 


o  w.  —  Du'qo';  —flic. 


-=— D(  SQa.  —vqS,  —V; )  —2  3  u\  —  Hu.  s. 


a  "5, 

— r- — rl^o,  =— A’Io’qU.  —  Ligo.  — —s.  U.  ). 
a  t 


B.  Special  Cases 

At  this  point  certain  special  cases  will  be  considered  which  can  be  si 
analytically. 


1.  Weak  Field  Case, 

For  this  case  the  atom  remains  near  the  ground  state  (w~— l)  and  we  can  ■ 
!r0=— 1.  w.  =0.  The  remaining  equations  are 


-=—3ut  -f  Ar,  —Vis, 


-=— Au.  —3  v. , 


— f — <  Kv0~ir)$.  =-Ki$0u, -v , ']. 

dr 


!  3.20c  - 


The  terms  on  the  right  in  equation  (3.1Qd)  are  linear  in  12.  [I\ 


.-n  •->, 


i\  =n  a ■  .m  ; ... 


Also.  since  u0  is  linear  in  H.  Ku0  is  quadratic  in  H.  The  dominant  term  in  equa¬ 
tion  (3.20c  i  is  irs..  A  solution  for  ■$>  of  the  form 

st  =Sj( t  )cosvt  -s,{t  )sin  ut  (3.21  i 

is  assumed.  If  coefficients  $\{t)  and  s-:(t)  vary  much  more  slowly  than  cos vt  and 
sin L>t  so  that  d2sL  dt'2  and  d2s.:  d t~  can  be  neglected  with  respect  to  dnds,  dt  and 
2 i/d s j  d t.  then  equation  (3.20c)  becomes 

d  s o  ds ! 

— cos ut - sim st)=—K  sQut  -~ut  -~u0(sxcosut  —s-sintd  V .  (3.221 

d '  dt 

Defining  if  =  ut—irt.  the  equation  for  £  is 


(3.23) 


The  solution  of  this  equation  is 


;(<)— n  I  ‘ 


st  Oo)^  to- 


(3.2-4) 


If  Sjidj  and  s.:ij)  vary  much  more  slowly  than  e  31 .  then  they  may  be  removed 
from  the  integral,  which  can  then  be  evaluated.  The  result  is 


sAn=-D,1n 


3cosut  ( 14-! )-!-( A-rldsinW  (1—0 
32- f(A-fi/)"' 


3cosnt  (1—i )— ( A— i/)sim /t  ( 1  —  i ) 

-v‘>  It  -  o 

d-+(A- i/j- 


— v2.s.,n 


dsini/'  f  1  — j )— ( A-ftOcosi/f  ( 1  -f  i ) 
d-.ftA-ri/i2  . 


«  *  >  ~  ^  a  *< 

'  -  V/V  fc  rt.ii  *>  w  ►  1 

kirk. 


■  '« , C.,’ ,i',v N 4'\ J 

*  '  *  h  hh  *.  k  A  fc  k  *  *  •  *  '  “  -  *  - 


3s\nut  i  1  —  i N— -  A— u'ccsut '  1—i  1 
3‘-A.  A—i/}~  J' 

Noting  that  Rebf’=y.  and  In'y=v,,  substituting  these  expressions  :'cr  u.  and  ", 
into  equation  (3.22).  and  equating  the  coefficients  of  sini/i  and  cos ut .  the  equation: 
for  and  are  found  to  be 

ds, 

- ——a  5  ,—75'.  1 3.26a 


- C o  S  y  , 


(3.25b 


where  a  and  7  are  given  by 


fl-Ar-n2 

3—{A—L/h0 

3 -A  A — idsq 

8my 

J^A-ur 

3--(A-»f 

a=AA_A - - — _A - A - -11  (3.27a; 

8m  u  i-q-{A— 1/}‘  3-—{A—u)~  \ 

_ fk2 fb’  d50— (A-v)  j.s0—  A-u)  ^  irk- Cl  o„v 

Smu  ,j--r{A— ic )-  d*— -A— uf  2m  u  0  (3._.  ^ 

The  solution  of  this  two-equation  system  has  the  form 

Si(t  )=e-a(  (.4  sin —Bccsut ),  ( 3.2Sa , 

5o(i  )=e~ai  (Bsiaut  —.4  cos ut ),  ( 3 . 2 S b  1 

where  .4  =— .s-:(0)  and  i^-s^O).  Thus  the  equation  for  s, ,  equation  (3.21).  becomes 

s(  (i  )=e~at  (.4  sin(7 —u)t  +B cos(';— u)t ' ,  (3. 25  ! 

which  is  valid  whenever:  l)  3>>a, 7  so  that  the  removal  of  s t  and  s.;  from  the 
integral  of  equation  (3.24)  is  valid:  2)  v>>a. 7  so  that  neglecting  d‘s.  d:2  and 
dT  d;:  in  the  derivation  of  equation  (3.22)  is  valid:  and  3)  3>>fi  so  that  t  ices 
not  differ  siznificantlv  from  -1. 


*i  '  l *  -1 


2.  Adiabatic  Approximation,  3  Is  The  Largest  Frequency 


A  more  general  case  which  does  not  place  as  strong  a  condition  on  F  and  t  is 
to  assume  that  3  is  much  larger  than  any  of  the  other  frequencies  in  the  problem. 
The  Bloch  variables  as  well  as  .5  are  assumed  to  have  the  form 

/  (  /  )=/0-r/  )cosut  ~f  2(t)s\nvt  (3.30 ! 

=/0+ReT(Oe,i"\ 

where 


When  the  second  derivative  of  5  as  well  as  the  nonlinear  terms  of  equa: 
are  negligible,  the  equations  of  motion  for  U ,  V.  It',  and  5  are 

^L=-(iu-~3)ir~AV-rC.s0W-rP.w0S. 
d  t 

—=-Ai'-{iu~3)V~Q\V. 

df 

^L=-nSoi' -nv  ~(iu^23)\v -nu0s. 


d  t  2u 


(SqL  4-1  4!i0O  ). 


(3.31 
ions  !  3.13 

i  3.32a  i 

(3.32b: 

(3.32c 

(3.32  i. 


We  consider  the  case  where  3  is  much  larger  than 


their  steady  state  values  much  faster  than  5,  and  at  any  time  the  values  rf  T.  f, 

and  H'  in  terms  of  5  can  be  approximated  bv  setting  bLIL  to  ger 

*  *  D  dt  dt  <i  t 


and  solving 


( !  u-~ 3)  — A  — fN'o 

’f] 

'  fiir.,  ' 

A  {iu—3)  —0 

V 

= 

0 

fh>0  n  1  i  u-r2  3) 

\\\ 

--“o 

When  these  values  for  U  and  \'  are  used  in  the  eauation  for  5  we  set 


dS 


d  t 


-(ar+i'7)S, 


wnere  ct  a 


are  given  bv 


Ot- 


fr/c2n- 


-4  mu 


Di.\\-D  o.V2 


D{+D? 


tk-a2 


Amu 


D  o.Vj-rZ?  j.Vo 


D?  +D; 


tk- n 


2  m  u 


and  .Vj,  A',,  and  £>■>  are  given  by 


A  i — 23~u.'qSq — 2dAu'Q-r-driug — z/~ 


ti’O-So, 


A  o — id i£ q  Air q+3 3 wqSq), 


Z)1=id5J2+A-+n2-^), 


D  o=2  3(  d-+ m) 


(3.: 


(3.3.1 


(3.3c 


(3.3( 

(3.3C 

(3.3i 

(3.3C 


As  can  be  seen  a  is  again  a  decay  constant  for  the  amplitude  of  o.  and  o  is  a  f 
quency  shift.  These  expressions  for  a  and  7  are  valid  whenever  they  are  much  1 
than  3  and  v. 


C.  Numerical  Results 


The  general  equations  (3.16)  derived  in  section  III. A  were  integrate]  mini*, 
cally  for  a  variety  of  initial  values  and  fr-quendos  using  a  fourth-or  Rr  Run 
Kutta  method  il  l).  This  allowed  comparison  to  the  anal'/tT-;.]  solutions  d-*ri\ 


here  for  special  cases.  The  Pascal  programs  used  for  this  are  listed  in  the  Appen- 


With  r=:/<  equations  (3.16)  take  the  form 


d  u  _  3  ^ 


a  n 

- V~, - IL'S  , 

V  V 


dr  u  v  u 


JT  SU  -  r  )  — — IT  —  1 

dr  u  ^  v  ^ 


1 3.37a  } 


3.37  b 


A'  ,  , 
{su+v)-s 
ir 


1 3.37  d 


(3.37e ' 


In  all  of  the  calculations  k  was  assigned  a  value  of  ICVitT1  and  m  was  assigned  a 
value  of  10-:°kg.  These  values  were  chosen  because  this  k  is  typical  for  visible 
light,  and  10-*skg  is  about  the  mass  cf  barium. 

Figure  12  is  a  graph  of  s  versus  time  for  i/=IOs.  i=10>.  A=— 105,  and  H— 5a104. 
Figure  13  is  a  graph  of  the  natural  log  of  Tie  amplitude.  A  =( 5TiU— .?„:n)  2.  of  s . 
versus  time  for  the  same  values  of  u.  3,  A.  and  C.  In  all  cases  the  Ln(A)  versus 
time  plots  were  linear,  indicating  exponential  decay.  A  decay  constant  given  by 


In  A  (/o)  —In  Altj) 


(3.3$) 


was  calculated. 

Equation  (3.27a)  predicts  a  value  for  a  when  f!  is  small.  Figures  11.  16.  and 
16  compare  a  calculated  from  equation  (3.27a)  and  a  obtained  from  numerical 
integration  of  equations  (3.37).  Figure  1!  is  a  plot  of  a  versus  A  for  ^=10".  3— 10‘. 
and  fl-- 5xl04.  We  see  that  the  cool  ini  rate  is  a  maximum  for  A:  —  r.  Abo  that 


■  O.  *JVa.TMW 


the  eccung  rate  oeccmes  negative,  wmch  corresponds  to  heating,  '.vhen  A  >0.  Fig 
ure  15  is  a  plot  of  a  versus  3  for  ^=105.  A=— 1.5x10°,  and  fl=5xl0'i.  Notice  tha 
maximum  ccciing  occurs  for  i~l2 u.  Figure  16  is  a  plot  of  a  versus  P.  for  i'=lQ- 
i=10°.  and  A=— 103.  The  conditions  under  which  equation  (3.27a)  is  valid  ar 
that  P  be  enough  smaller  than  3  so  that  w  does  not  differ  appreciably  from  — 1 
For  the  values  of  u,  A.  and  3  used  for  this  plot,  there  is  good  3.cr  m  p  jV; 

n<o,4i. 

Equation  (3.35a)  gives  a  value  for  a  when  3. u>  ><a,q.  Figures  17.  15.  T 
compare  a  calculated  from  equation  (3.35a)  and  a  obtained  from  numerics 
integration  of  equations  (3.37).  Figure  17  is  a  plot  of  a  versus  A  for  v=103,  J=10: 
and  n=l03.  The  cooling  rate  again  becomes  negative  for  positive  A.  Figure  IS  is 
plot  of  a  versus  3  for  ^=105.  A=-103,  and  H=103.  Again  maximum  cooling  occur 
for  3^v2U.  Figure  19  is  a  plot  of  <n  versus  P  for  z/=103,  i=10'.  and  A=— 103.  Fo 
larger  P.  the  two  curves  again  fail  to  agree,  but  the  value  for  a  calculated  from 


equation 


does  reach  a  maximum  and  begin  to  decrease.  This  se: 


behavior  agrees  with  that  shown  by  the  numerical  integration.  Equation  (3.27a 
does  not  exhibit  this. 


IW.V 


o.  versus  (3 


f  isurf  15.  Plot  of  a  versus  3  for  i>—  K>r’,  Ji- 
1  i n o  rfprr-f'nts  flu*  vnlnr-?  of  n  onloulnt r<l  from  i 
v:iiiifs  of  o  raP-uhitfil  from  equation  (5.571. 
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IV.  Single  Atom  Spectroscopy:  Ladder-Configuration 


1 

r 


m 


For  the  four-level  svstem  shown  in  Figure  2.  the  following  assumotions  about 


the  Einstein  A  coefficients  and  the  dipole  transition  moments  between  the  various 


states  will  be  made: 


-4  30  -d  31  d  ;q=0. 


-d  3 1  <  <  -d  10.-d3.3- 


(-4.1b) 


which  imply 


^30-^31— 


(4.2a  1 


/A  1  <  <  MlO'^32- 


(4.2b! 


In  addition.  (il0,  n-2v  and  1*3-2  will  be  taken  to  be  real.  The  applied  electric  field  is 


assumed  to  have  the  form 


E(t)=E icosojf  —  FTcos-'of  — Eicos-gjf 


with  components  at  frequencies  ~v>,  and  ~3  near  to  resonance  with  the  th 


transition  frequencies  w10.  w;i.  and  w3;  respectively. 


A.  Dressed  Atom  Approach 


A  simple  picture  which  shows  some  of  the  qualitative  aspects  of  th”  la 


configuration  is  the  dressed  atom  approach.  Without  spontaneous  emis-ru'..  tin. 


Hamiltonian  matrix  for  this  svstem  in  the  rotating  wave  at* proximal!; 


diagonal  (21)  when  the  basis  states  are:  1)  the  atom  in  the  ground  state  an-!  •;  • 


photons  cf  frequency  Wj,  n_.  photons  of  frequency  c,.  and  n6  photon.- 


•rw  r*  r 


-.'3,  state  |0.rc..no.n3>:  2)  the  atom  in  the  first  excited  state  and  n.— 1  pnetc 
frequency  etc.,  state  1 1 . n  j  —  l.n.:,,n3>:  3)  state  |2, n  t— 1  .n.2— 1  ,n3> :  and  -! 

[ 3. n t  —  l.n.:— l.n3— 1>.  A  representative  block  is 


Hlj=n'nl~l+n.z~.:-ni~2i 


0  V2i'Ol  0  0 

-fcjfij  -Aj  hr  a  o 

0  -VziQ,:  —A,— A,  Vi  itt J 

0  0  -hi  Qj  -Aj-A2-A3 


where  /  is  the  identitv  matrix  and 


Aj— U.-J— WjQ, 


-Xo - (jjfs - 


2  ^21  > 


■^3  "3  "32 


are  the  detuning  n  the  various  transitions,  and 


VioE  i 


^3 1  2 


M  32^3 


are  the  on-resonance  Rabi  frequencies  for  the  three  transitions  being  consk. 


This  is  a  nearly  degenerate  system:  the  differences  in  energy  between  the  5 
being  just  A  times  the  appropriate  A,'s.  The  dressed  states  are  the  ; 


v£y  «lk  ■C-? -If".'’  v 


obtained  bv  diagonalizing  the  Hamiltonian.  The  dressed  state  enerait 


E  '  =E  —~a\ 

•‘-I  '  ■*  > 


where  the  X,-  are  solutions  of 

X(X4-A1XX-A1-A:XX-A14-A,-A3)-lr>.fX(X-A1^A:-A3) 

4 

4  4  lb 

For  f>:=0  and  A1;  A3=0  the  solutions  to  equation  (4.S)  are 


X^-HQj  (4 

X2=+V2fX  (4 

x3=-y2o3  (4 

x4=+^n3  (4 


These  X’s  correspond  to  eigenstates  for  the  characteristic  block  of  the  Harniltc: 
matrix.  The  transformation  matrix  between  the  old  states  and  these  eigenstate 


!  1  0  0 
1)00 
0  0  !  1 
0  0  1  i 


(4 


This  transformation  diagonalizes  the  Hamiltonain  when  Aj.  A3— 0  and  fT=0.  \Y 
the  Hamiltonian  is  diagonal,  the  basis  states  are  stationary,  that  is  there  a;\ 
transitions  between  them.  Figure  20  is  an  energy  level  diagram  for  the 
states.  If  Ho  is  turned  on,  there  will  be  transitions  between  these  s' at. vs.  In  •  hm 


ystem.  the  transition  frequency  between  state  1 1  ,n  (  —  1  .n, :  -  and  |2.  \ ,  —  1 .  r3  ■ 


-Uj.  There  are  four  transition  frequencies  among  the  new  stat 


TiH,— n.).  Thus  th>’  dressed  state  approach  pr-vid- 


bv  AHiTh-fb)  and 


A 

y. 


«p 

>* 

«.■ 


dec  ail  what  the  cha.ra.cter  of  the  hucresoer.ee  will  be.  it  is  heioiui  to  turn  to  tr.e 
optical  Bloch  equations. 

B.  Optical  Bloch  Equation  Approach 

The  strategy  we  will  use  to  extract  numerical  information  about  the  flucres- 
cence.  such  as  the  average  dwell  time  in  the  upper  two  states,  is  similar  to  the 
approach  taken  in  reference  (23).  Reference  (23)  deals  with  coherent  excitation 
the  V  configuration.  In  that  paper  the  authors  show  that  under  certain  general 
conditions,  excitation  of  the  weak  transition  can  be  treated  as  a  rate  process  when 
viewed  over  coarse-grained  intervals  of  time.  Being  able  to  treat  the  excitation  . 
the  weak  transition  as  a  rate  process  greatly  simplifies  the  analysis  of  the  statistics 
of  the  fluorescence.  A  similar  proof  will  be  given  for  the  Ladder  configuration. 


1.  Development  of  Optical  Bloch  Equations 


The  equations  of  motion  for  the  components  of  the  density-matrix.  p.,„.  of  a 
coherently  excited,  multi-level  system,  in  the  electric  dipole  approximation  are  (27 


d  t  dt 


=  -  Tj  AnrPn+H  Pr  +  ~T~TjiPnr  Urn  Prn  > 


r  <  n 


r  >  n 


itr 


for  n  =  m.  and 


d  P nrr. 


d  t 


-Va  V  -l  ,.i:,  V  4  ' p  -uibliiv' ,  a  _ 

* _ i  '  »;r  /  j  -  1  *nr  t  n vi  >  >\r  r*- 

HI 


for  n  r-  rn.  where  Pn  is  the  probability  to  1-  in  'ho  :rh  -*•> 
elements  or  coherences  varv  much  more  ra;  id’.v  than  or 


i-l.lla! 


i-l.l  lb 


iTiTt 


venien:  :o  work  with  the  slowly  varying  coherences.  7,,.  denned  by  th 


Poi—^qi? 


1 4.12a ) 


pQ-2 — aZ2e 


(4.121 


POZ — & 03 e 


>  ( i  +w3-h«’  3 ;  t 


(4.12c) 


Pi*  j  o  € 


(4. 12d) 


P 13— <713e 


i4.12e! 


p.i3 — a  .^e 


(4.12f) 


We  shall  also  work  with  the  inversions 


u’i  —P  i~P  o> 


(4.13a) 


W2=P  2~P  l> 


(4.13b) 


WZ~P  Z~P  2; 


(4.13c) 


and  the  probabilities 


P.=P0+P  i, 


(4.14a) 


P+=P2+P3. 


(4.14b) 


Note  that  P _  is  the  probability  to  be  in  one  of  the  lower  two  states  and  P_  is  the 
probability  to  be  in  one  of  the  upper  two  states.  Consider  two  cases:  l)  when 
P_— 1  and  F\=  0  we  would  expect  to  observe  scattered  light  of  frequency  a,.  2' 
when  P _= 0  and  P 1  we  expect  to  observe  scattered  light  of  frequency  -g.  During 
an  actual  experiment,  the  status  of  the  system  will  be  constantly  monitor--!  a 
time  scale  which  is  short  compared  to  the  rate  of  change  cf  /’_  am!  /’_  -1-rive- i 
from  the  Schrodinger  equation,  so  that  case  l)  or  2)  will  always  apply,  and  th- 
system  will  change  abruptly  between  the  two.  This  is  what  is  tm-a:;t  by 


3<713 
d  t 


8 


_  i„  j  __i  \  _  \  ,t  o  -  -  _:,,o 

2!  --i  iq  — jo  l  .  — j  — 1;  13  -  \  -  Zi  —  2  1  -  -ju  j3  —  -  ■  ■  ■ 

dcr.,-, 

=—  +i  .A  32 — .4  ,,  !  — :  A,  <7.y,~ '-2J  — '-2i  B->  r 


d£ 

dP_  dP 


d£  d£ 


-=’"2.4  .>j(PT— >rz)-~v2i  B_o (Toj — ri: 


4.15 


i .  1 5 : 


^  5 

.i  *  - 
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I* 


r 

* 


■r. 


i. 

fc  5 


2.  Steady  State  Solutions  of  the  Bloch  Equations 

We  will  first  consider  the  steadv  state  solutions  to  equations  (4.15).  For  a 


braic  simplicity  the  detunings  on  the  strong  transitions,  \  and  A3.  wi 


will  be  set 

zero.  We  will  define  the  complex  quantities.  A'x.  Ah.  .V3,  and  A’4  as 

A  j=.4. 3i-f2:  A:,  (4.16a 

A  2  =.4  32-1-2 :  A3.  (4.16b  ; 

A3=,4  iod-,4 2i“r2i  A3.  (4.16c) 

N^=A  10+-4  32+2 iAj.  (4 . 1 6d  ; 

The  system  that  must  be  solved  for  the  steady  state  solutions  of  the  cq-'s  in  terms 
of  the  u’i’s  is 
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itr.U' 
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(4.17 
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v.  here  D .  the  determinant  of  the  matrix  in  equation  .4.17'.  is 
D  =rJ.Y,.V3-PJU  ,0;  .V1.V,.V3^oi:A-^o2.\-3, 

^r>|.4  34.v;V3.v4-o^v3_o^A-:)_o^  ;0.43o(n;--n|-.v,.v4; 

-n32.4 10.4  32(^-^f -^-v3-\'4)^  10-a  32.v,(  .v;.y3a>o^y3-q|.y  , . 

In  matrix  form  equations  (4. IS)  are 


*01  rll  -12  -13  ’x'l 

*12  =  -21  -22  -23  '-f  2 

*23  -31  -32  -33  ’x3 


where  the  q,  's  are  defined  by  equations  (4. IS).  Taking  y. =Im(  r. ,  :.  one  mav  writ- 
the  equations  for  the  steady  state  values  of  the  w,'s  as 


1  ~  !  o — 2  n !  y  x  J  -r  Ho  y  .1 J  )  U’  J  U-(  .4  !  o — 2  n  J  y !  O — n2 !/ . V,  j  U1  .J 

'r(^-d  io--^iy  13d-n.2y.23)  w3=^(A  iQ  —  id-d  21) 


(4.21a) 


(*d  10  3.4  21  ^33)  ■  y  1 1  2noy-j1-;- n.jy31  uq 


-M—A  io~ 2.42l-f.4  3 2 ) d" H j y  r: — 2 flj y o — fl 3 y3 2 ’  u ■ 


>  A  OI  k  \ 
t  4  .  _  1  0  f 


~^(  -d  10d-2.d  21-3.4  32)+^^  !3  2fl.2y23-rfl.3y3 


33  ^  3“'  2  —  ’ili.-l  in— .-1 


3'  .  •  '*  ■-  .%  ,*•  ■,  \  •«  •»  %  *.  %  *.A  A  A  Jj 


or  in  matrix  norm 


all  3  12  3 13  u'i  -d  io— ’“*-4  21 

3  21  3 22  3  23  ’r  2  =  ^  -d  21  —  ^(-d  10_r-4  32) 
3  31  3  32  3  33  ’3'3  -^32— ^-^21 


where  the  a,-, 's  are  defined  by  equations  (4.21).  We  define  .4  as  the  dote 
of  the  a,j  matrix.  The  solutions  for  the  w,  's  are 
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Pascal  program  listed  :n  the  Appendix.  FA  '.vas  calculated  as  a  Auction  of  A: 
various  values  cf  fA  W,  and  123.  Figure  21  is  a  plot  c-f  FA  versus  A._>  for  ,4;3=10’ 
A  ji=1  .  A  3o=l  .25x10°.  ^=10'.  1A=103.  and  123=2x10'.  Note  the  peaks  in  the  spec 
trum  cf  FA  at  A:=— j — 1A)  and  A:=— ’-21 12j— O3)  which  agree  with  the  dresse 
states  results. 


3.  Adiabatic  Solution  of  the  Bloch  Equations 

A  more  useful  approximation  than  the  steady  state  solution  is  the  adiabat: 
approximation.  To  determine  the  statistics  of  fluorescence,  we  are  interested  i: 
how  FA  and  P_  vary,  not  so  much  in  how  the  coherences  and  inversions  on  :h 
individual  transitions  vary.  Note  that  in  the  equations  for  the  time  derivatives  c 
u‘-r  w3,  <701.  <703.  <712,  <713,  and  <j.,3  there  are  large  self  relaxation  terms.  The  proba 
bilities  FA  and  P  +  vary  much  more  slowly  than  these  other  quantities.  At  t hi 
point  it  is  difficult  to  tell  how  a02  varies.  It  is  reasonable  to  suppose  that  thes 
quantities  relax  to  their  steady  state  values  much  faster  than  P _  and  P _  change 
We  begin  by  solving  for  the  steady  state  solutions  of  <t01.  <t12.  and  a in  terms  o 
’iq.  w3.  T=7:o— -7 Q.:.  r=i(-7;0--0.:).  and  IV=P,-P_.  The  equation  for  <r01  can  b> 
solved  by  inspection.  The  solution  is 
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The  adiabatic  solutions  for  <jv2,  cr.^  and  cr03  are 
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where  £>,.  t  he  d'-tenninant  of  the  matrix  in  equation  (4.20; 
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equations  ;  or  t r.ese  quantities  are 
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Putting  the  steaiv  state  solutions  of  »r.  and  w,  into  the  eouatior.s  her  the 


and  putting  the  expressions  for  the  appropriate  parts  of  the  at.'s  into  the  equa¬ 


tions  for  U .  P.  and  H'.  ws 
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where  the  constants  are  given  bv 
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I.  Conclusion 


A.  Discussion  of  Results 

In  summary,  the  problem  of  laser  coding  of  a  single  ion  or  atom  in  a 


monic  trao  was  considered.  A  simple  rheorv  of  sideband 


developed.  In  the  limit  that  the  particle's  secular  motion  car.  be  treats :i  send 
sically,  the  theory  allows  the  calculation  of  a  logarithmic  cooling  rate  e 
numerically  or.  for  two  special  cases,  analytically. 

This  theory  of  laser  cooling  differs  from  previous  theories  in  that  this  t; 
predicts  cooling  rates,  as  opposed  to  minimum  temperatures  of  cooling,  f 
advantages  of  this  development  are:  li  It  provides  a  simple  physical  picture  : 
cooling  process.  2)  No  assumptions  needed  to  be  made  in  the  derivation  of  * 
tions  (3.16).  the  equations  at  the  core  of  theory,  about  the  relative  mu 
the  radiative  decay  constant  3,  the  oscillation  frequency  in  the  trap  v.  the  d> 
ing  A.  or  the  Rabi  frequency  Cl.  That  is  the  theory  as  whole  is  not  a  weak 
■:  fid  <  3'  theory  or  a  strong  trap  (u>  >3)  theory,  etc.  3)  Calculating  coding 
for  a  set  of  parameters,  even  in  a  region  where  neither  of  the  analytical  sdu 
applies  is  a  relatively  simple  numerical  procedure. 

As  has  been  mentioned,  other  theories  con  cent  rate  on  deter::::::::;.: 
tetnj  era*  nr- s  f  coding,  net  coding  rat  vs.  It  a  no  and  Wineiani  :  1  s  :  .  j  _  ; 
use  in--. her-  at  mt-mgy  rate  equal i -ms.  among  ot  imr  nvdt-.-ds.  ’  ■  r 


minimum  ■•■mperat  ur-v.  They  make  tin-  weak  fi-M  approxiu 
is  iti'erest  ;:m  P >  c  . : :  o r--  f ■ : u  ct i - ■  r. o !  f  r;:i  Coir 


►  n 


equation  !  10ai  of  reference  US),  to  equation  (3.27a:.  the  equation  for  a  f 
weak  held  case.  Xote  that  although  a  is  the  decay  constant  for  the  amplif 
the  motion,  the  energy  of  a  harmonic  oscillator  is  proportional  to  the  a:r.: 
squared,  so  that  the  decay  constant  for  the  energy  should  be  twice  a.  Their 
tion  is 


— 2(*b-oWi'-<  v,-> 


o~'- 


v— ■ 


In  the  notation  of  reference  (IS).  7  is  the  natural  linewidth.  m0-w  is 
R  =1 lik  f'  2 M.  fi3  is  the  scattering  coefficient,  equal  to  1  3  for  isotropic  scatte: 
and  7,.  the  average  scattering  rate,  is  given  by 

7,=^r-<T~'w)>, 

where  a\  :  w.T)  is  the  Doppler-shifted  natural  line  shape.  /  is  the  intensity  of 
cooling  laser  which  is  proportional  to  the  square  of  the  Rabi  frequency.  fh 

The  first  term  in  equation  (5.1)  describes  cooling,  while  the  second  t 
73  R  ( 1  -r/?: )  corresponds  to  recoil  heating.  The  energy  of  a  harmonic  osciilat; 
proportional  to  the  average  over  one  cycle  of  the  v:’  city  squared,  so  that 
cooling  term  is  linear  in  the  energy,  and  corresponds  to  exponential  ccciin 
predicted  by  the  present  theory.  The  resonant  denominator  is  the  same  for  ! 
the  weak  field  expression  for  a  and  the  cooling  term  in  equation  (5.1).  Also, 
cooling  term  in  equation  (5.1)  scales  as  fh  as  does  a  for  small  fh  The  depend 
of  a  on  3  and  u  differs  from  the  cooling  term  presented  by  It  a  no  and  Wine! 
In  choosing  o  from  equation  (3.27a)  or  the  expression  presented  by  \ Vineland 
llano,  we  expect  the  ex;  ression  for  a  to  be  the  appropriate  value  f  :  an  en 


decay  cc 
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laser  cooling  because 
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The  preset: 


'laser  intensity’1;.  D:  and  to  some  extent,  by  cheesing  appropriate  transitions.  the 
radiative  decay  constant.  3.  An  appropriate  transition  for  ccciina  tv :.uid 
between  the  ground  state  and  some  excited  state,  with  no  possible  decay  rath 
back  to  the  ground  state  except  straight  down,  so  that  potential  cooling  tin:-.-  is 
net  wasted  with  the  atom  sitting  in  intermediate  levels.  Such  an  ideal  situation 
can  be  hard  to  cone  by  in  practice.  The  frequency  of  the  cooling  transition  must 
also  be  accessible  to  a  tunable  laser.  An  additional  consideration  in  selecting  a 
transition  is  that  the  rate  of  cooling  transitions  per  second  is  limited  by  3.  so  :ho‘ 
it  would  seem  that  large  3  would  be  better.  However,  when  3  becomes  Inroe 
enough  that  the  sidebands  are  no  longer  well  separated  from  the  central  teak 
cooling  is  slowed  by  transitions  between  the  sidebands.  The  maximum  cooling 
rate  occurs  for  3^V2U.  Once  an  appropriate  transition  has  been  chosen,  v  should 
be  adjusted  as  much  as  possible  so  that  This  is  not  always  possible:  in  the 

experiment  performed  by  Xeuhauser,  et  al  (31),  3  was  2 tx  19MHz  and  v  was 
2rrx2.4MIIz.  For  well  separated  sidebands,  the  value  of  A  for  the  maximum  ced¬ 
ing  rate  is  just  —v.  If  the  sidebands  are  not  well  separated,  the  value  of  A 
corresponding  to  the  maximum  cooling  rate  is  shifted  toward  zero.  For  giver, 
values  of  3  and  u.  plots  of  a  versus  A  are  easily  generated,  even  in  redoes  where 
neither  analytical  solution  applies,  and  the  best  value  of  A  can  be  taken  from  a 


graph.  The  same  is  true  of  fl. 

The  spectroscopy  of  a  single  atom  in  the  ladder  condign  rat 
treated  theoretically.  A  dressed  atom  approach  was  used  to  provi 
information  about  the  system.  It  was  shown  that  the  weak  transude 
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Appendix:  Programs  Used  for  Numerical  Calculations 


I.  Numerical  Integration  of  Equations  3.37 

Program  PNEW2; 

{*>} 

{SI  GRAPH.?} 

LABEL  1.2: 

CONST 

N:integer=  100: 

hbar=1.055E-3-l: 

mass=lE-25: 

k=lE-07: 

nu=lE— 06: 

beta=10: 

delta=-l: 

{  omega— 1:  } 
dt— 0.05; 


o 

r 


A  /•. 

'X‘  c 

vN  r- 


tj 


TYPE 

BUFF=ARRAY[0..406]  OF  BYTE: 
number= array  [1.. 250]  of  real; 
gob=array[l..250]  of  integer: 
rungekutta=array['u'..Vj  of  real: 
title=string[lo:; 

•—’l  j 

MAR 

BUFFER:BLTF: 

SAMPLE  :BYTE: 

U.LOL.AdXTEGER: 

DATA:text: 

All  AMA  W.AX.C.X  l:  n  u  m  her: 

U.VAMX.Y, KAPPA,  omega:  real: 

K 1 ,  K2 ,  K3 ,  K-l :  r  u  n  g  e  k  u  1 1  a : 

PROCEDURE  R  U  NG  E 1 1 1 .  v ,  \v .  :<  .y .  k  a  p  p  a :  re  a ! :  MAR  rk:rungekut  ta ’ 
BEGIN 


»« 4***1 


£ 


I 

? 

I 


S' 


v, 


rk.’u' 

rk;'v" 

rk.'w' 

rki’x" 

pk  y 

END: 


it  ‘f-bet a*': — ieita'v — jmeca*,.vtx  : 
<it*i  — iei  t  a  *  a- ct  a  *  v — c mega  ‘  \v  : 

:dl  *  '.-c  m.e-ta  *  x*u — v  -2  *  bet  a  *  •  ’.v  —  1 '  : 

-i .  *.  i 

dt*i -kappa*  x‘ x-v  -x  : 


PROCEDURE  \  1 1  NIT  IU\  [  l  n  u  m. :  a  u  m  b  e  r : 
BEGIN 
min:=num  ly. 

FOR  I:=2  to  L  DO 


Liinteat 


BEGIN 

IF  num.I’  <  min  THEN  min:=numT: 
END: 

END: 


\'AR  rnimrea)  : 


PROCEDURE  M-YXIMlNIinummuruber:  L:i 
BEGIN 

max:=num;  l'; 

FOR  I:=2  to  L  DO 
BEGIN 

IF  num[I]  >  max  THEN  inax:=rm:nT: 
END: 

END: 


n  tetter: 


WAR  m 


PROCEDURE  SCALE(min.max:reaI:  numtnumber:  UUttger:  \’AR  mark:g : i  : 
V.AR  s:real; 

BEGIN 

s:=(max-min)/lSO: 

F OR  I:=l  to  L  DO 
BEGIN 

m  a  r  k  [  I  ] : = r  o  u  n  d  ( ( n  u  m  [  ]]■ -  m :  n )  /  ~ ) : 

END; 

END; 

PROCEDURE  PRINT( nitm :n u m ber :  L?flas:inte-zer:  Tx.itle:  \'AR  :  U  : 

WAR  min.maxmeal:  mnrk:gob: 

BEGIN 

GRAPIIMODE: 

D  RA\  V(  10. 100.200,100. 3): 

DRAW(  10.100.10.10,3); 

D  RA\V(2  60. 1 00 .2  GO .  1 0 . 3 ) : 

DRAW( TO. 10.200.10.3): 


,s2 


FOR  I:=  1  TO  9  DO 
BEGIN 
J;  =  1G0-I‘1S: 
DRAWjlO.J.lS.JM: 
END: 


FOR  I:=l  TO  9  DO 
BEGIN 
J:  =  10— 1*25: 
DRAW(J,1G0.J.1$7.3;: 
END: 


MI  NIMUM(  aura.  L.  rain  ): 
M\XLMTDvI(num.L.max'): 


IF  flag=S  THEN  par:— 0: 

IF  flag=2  THEN  p ar: = ( num 'L j  - nu m !  1  j )* n u/ ( dt * 3 * ( N *  1 E  — 02 — 2 5i 
IF  flag=l  THEN  par:=-LN((max-min)/2): 

GOTOXY(1.25): 

WRITE(min:lO.'  \max:10,'  ’.(max-min):10); 

GOTOXY(IJ): 

\\EJTE(((max— rain'I/2):10.'  \T): 

IF  flag=2  THEN  WRITEf  ’.par); 

SCALE(  min.  max.  num.L.  mark): 

PLOT(ll,l90-mark[l],3): 

EF  flag>l  THEN  CERCLE(ll.l90-mark;l’.3.3); 
J:=ROUND(250/(L-1)); 


FOR  I:=l  TO  (L-l)  DO 
BEGIN 

PLOT(llE(lNJ).190-mark[I-l:,3): 

IF  flag>  1  THEN  CIRCLE(11P!  W)JG0-mark[lEl\3.3): 

IF  flag <3  THEN  DRA\V(ll-i-(I-l)*J.190-mark;i;.ll^I‘J.100-raark;i-i;.3) 
END: 


IF  flag>l  THEN 
BEGIN 


FOR  I:=l  to  SO  DO 
BEGIN 

GETPIC(BlTFER.i*(I-J  ).0,l*il-l  }.200 ): 

\VRITE(LST.#27/@\#27,T): 

\VRITE(LST.#27,#7o.#20I.#0): 

FOR  .J:=  1  to  201  DO 
BEGIN 


V 


SAMPLE  :==3LTFER  J-i  -M 
WRITE.  LST.CHE  SAMPLE 
END: 

END: 

\\"RITELNiLST.==l l.=l l.=l  1  .=10. 10. =27. 'M '  ; 
END: 

ENT): 

FL  NOTION  KUTTAi rkl.rk2.rk3.rk-l:real):rea!: 
BEGIN 

KLTTA:=(  1  M  * ;  rkl-2  *  rk2-2  *  rk3-rk-l ': 

ENT): 

BEGIN 

FOR  LOL\:=l  to  5  DO 
BEGIN 

omega:=LOLATO: 

l':=0: 

Y:=0: 

\V:=-1 : 

X:=0.2: 

V:=0.2: 

K.APPA:=hbarJ'sqr(k)*cmega/(nu*2*mass); 
writelni  lst.’nu=  \nu): 
tvritelnf  1st. ‘bet  a=  \beta‘nu): 

\v  ri  t  el  n  ( 1st  F  del  t  a=  ’ .  •  .1  e  1 1  a  *  n  u } : 
v  ri  t  el  n  ( 1st . ’ o meg  a  =  '  .omega*  nu): 
w  r  i  t  e  1  n  ( 1  s  t .  ‘  d  l  a  u =  ' .  d  t  j ; 
w  ri  t  e!  i.  ( Is  t .  N .  ’  x  1 00  ’ ; : 
w  ri  t  el  :t :  1st. ==10.  A]  Ox 


FOR  1:=1  to  100  DO 
BEGIN 

RUNGE(L'.\'.W.X.Y.KAPPA.K1): 

RUNGE(I’-r-0.o*Kl  u’  A'-r-O.o  ,KL'v';.\Va0.5*K1i'v.-‘,.X— O.o'Kl  X  . 
V— 0.o*K  1  [V  j  Jv.-VPPA.K2): 

RUNGE(l'-K).o*K2;V;  .Y--0.o  *K2M\  .W-0.5  *K2''w'‘.X-  O.r,  •  K2  . 

Y-f0.5*K2[y;  .NAPPA.K3); 

RUNGE(U-FK3['u‘{,VAK3[’v'].\V-rK3['w'j.X-bK3i*'x'jA'— K:-)i’-y‘;.KA[’I'A. 
U : = U -f K U T  T  A  f  K 1  ['u'].K2  'u\K3|’uNK 1 N'; ): 

v:==\>iaTTA(Ki;N--;.K2:l’v';.K3i'v,;.ivi[,v,;}: 

W:— \V -f-KUT T  A  [  K 1  !'  w ' 1 .  K2  " u  " . K3  ’ Ev ‘ L K  ! ' '  M "  ' : 


X:=\— KL  TTAi Kl  x'  ,K2  ,K3  'x'.Kl  ’x"  ; 

Y:=T  —  KL'TTA, Ki  'y".K2  'y".K3  V  ,K t  : 
END:  {End  of  Outside  Lcoc} 


FOR  I:=l  to  250  DO 
BEGIN 

RUNGE(U.VAV.X.Y.KAPPA.Kl): 
RUNGE(U-y0.5‘Kr'u';.\'-0.5*Kl;'v';.\V-0.5*Kl'dv'  .X- 
Y— 0.5*Ivl  [V  j.KAPPA.K2): 

RUNGE(U-r-0.5‘K2:'u';A'-0.5'K2;'v'yW-0.5'K2  Nv".X- 
Y-0.5*K2[y;.KAPPA.X3): 

RL'NGE(L'— K3[Vj.Y— K3[  VDV— KSXvXX— K3  {x’.Y-i\3 
E':=U-rIvUTTA(Kl.'u’;  J\2:’u’].K3[‘u‘l.K-r‘u'i): 
v:=\  ainxtta(ki;'v';.k2;-v-;.k3:'vXk-B'v-;'): 

W:=\\  — KL'TTA! KrNv’;.K2['w’].K3[Nv’!.K-4’Nv"  ■  ; 

X: =X-r  KL'TTA ( I\  1  ■  ’  x  ’  LK2 V : .  K3 : '  x ' ; .  K-l ; V ' :  ’ 
Y:=Y-IvUTTA(Klf-v1.K2?VdK3l'vXK4''v'V;: 

AU[I;:=U: 

AV;ij:=V; 

A\V:i!:=\V: 

.\X;ih=X: 

END:  (End  of  Outside  Loop} 

PRINT!  AX. 250,  l.'Sl'.CTP: 


FOR  I:=l  to  N  DO 
BEGIN 

FOR  J:= 1  to  100  DO 
BEGIN 

RUNGE(UA'.\V.X.Y.IX\PPA.K1): 

RUNGEi'U— 0.5*Kl[’u’j,Y-F0.5cKl['v’].W— 0.5*Kr,-.v".X-o,5*Kl  "x". 

Y-f  0.5  *K1  [V  j  ,KAPP  A.K2); 

RLNGE (L  -r0.5*I\2[  u  ].\  A0.5*K2[v'y\V-y0.5*K2  ’\v':.X— 0.5*K2  x  . 
YA0.5*K2[y].KAPPA.K3): 

RUNGE(U  aK3  [‘ u PK3  [!  v'l . \VaK3  I' v.- '  1. X—K3 O’— K3 \K  A:  M  'A  ,!\ 
L  :=L  —KL  TTA!KF'u'yK2i'u':.K3  u'  .K-l  'u'  y 
Y:=V~Ki;TTAiKY'v‘!.K2:'v'!.K3''‘v'j.K4S-d*): 
\V:=\VaKUTTA(Ki;Nv';.K2Xv3.K3;N.';.K!Xv''  : 

X:=X-KUTTA(K1  S:-;.K2:NNK3rv';.K  F'v{  : 

5  :=A  A K U T T A ( K 1 1  ’ y ' y K 2 i ' y ' i . K3 , ' y ‘ j . K ■  1 : ' y  ‘  y 
END;  {End  of  Inside  Loop} 

END:  (End  of  Outside  EooiP 


vs 


* 


r  OR  I:=  L  to  250  DO 
BEGIN 

RUNGEI  LLYAV.XXLKAPPA.Kl): 

RL'NGEiL*— 0.5‘Kl  ALLY— 0.5 ‘Kl  "v  "  AV-0.5  ‘Kl  Nv'LX-0.5  *  Kl 
Y~ 0.5‘Kl  ,'y  XKAPPA.K2 

RL  NGEi  l  — 0 . 5  ‘ K2 !  V , . Y— 0 . 5  ‘ K2  L v ' '  A V— 0. 5  ‘  K2  Lv '  LX — 0 . 5  *  IX 
Y-0.5‘K2  [y;.  KXPPA.K3): 

RUXGE^U— K3[’u‘!,Y— K3:‘v'lAV— K3i”’.v':.X— KolAp’.Y— K3’'v".I' 
U  :=L’— KLTTAl’Kl  [  V]  ,K2 ; '  u’LK3  ['  u’*j  .Iv4  [V  j: 

V:==Y— KL'TTAlKl  •  V]  .K2|  V  j  .K3i  V;  .Ivi  [VI ): 

W  :=\\  —  KL  TTAiKl i’tv'j.IX  ' *,v ' j . K3 • ' \v ' I . K-4 ; '  V  ’  i : 

X:=X— KL'TTAlKl  :  V.J\2[Y;iv3[Vj,IvLV:): 

Y:= Y— IXUTTAiKl :  VLK2;  Vj,K3i  VLK-i  ’V'): 

Au;n=u: 

A\’il::=V: 

A\v:i;;=\v: 

AX:L:=X: 

EXD:  {End  of  Outside  Loop} 

PRINT(AX.250.1.'S2LC[2j): 

FOR  I:=  1  to  X  DO 
BEGIN 

FOR  J:=l  to  100  DO 
BEGIN 

RUNGEi  LLY.W.X.Y.KAPPA.KIj: 

RL  NGE(L  —0.5  ‘Kl  {u]  A— 0.5‘Kl  i  Vi  AY— 0.5  *Kl  [VLX— 0.5  ‘Kl 
Y-0.5*K1  i'y-:. IG-VPPA.IX): 

RLLXGEfL'— 0.5  *K2  ['u’l.Y-f  0.5  ‘IXL'v']  AV-r0.5*K2[‘wL.X— 0.5  *N2 
YA0.5‘KNyi.KAPPA.K3); 

RL’NGE(U— K3'-u']A'— K3[‘v'i  A\VK3[’'.v’’.X— K3’'x’i.Y— K3’-yYK 
LL— L'~ KU TTA(Kl [’u{.K2  'uAK3N.i'’.K-l,'u{): 
V:=\'-hKUTTA(Kl['v’i.K2i'v';.K3[‘v’j.K-i.’v'!  j: 

W  :=\\  +KL  TTAfKl  .K2[*\v”  j  .K3;  V  J  .K-l .  tv' 

X:=X-fKTTTA(Kl  [  Y).K2[  Y]  .K3rY]  .KEY'  '): 
Y:=Y-KLTTA(KI  I'y'i.K2{y',|.K3j'y'’.K3 
END:  {End  of  Inside  Loop}1 
END;  {End  of  Outside  Loop} 

FOR  L=1  to  250  DO 
BEGIN 

RUNGEI  L'A'AV.X .  Y.KAPPA .K 1 ): 

RUNGEfU-0.5‘Ki;'ir;.Y-0.5-KU'v';.\V-0.5'Kl"..  ,\  . 

Y—0.5 ‘Kl  ’  VLKNPPA.K2  y 


I  i 


.  -X -X  »*,  -X  -ov-X 

-•  ■A*.  if .-A. 


^V»  ■*-  ■-  -- »-  — 


RUNGE(U+0.5*K2[’u,],V+0.o*K2[V],W-r0.o*K2[,w,].X-r0.o*K2r'x,!. 

Y+0.5*K2[,yJ],KAPPA.K3); 

RUNGE(U+K3[’u,jlY+K3[V],\V+K3[,w']>X-rK3[’x'].Y-rK3[,y-’.K.\PPA.K4'.: 

U:=U-rKUTTA(Kl['u'],K2['u'].K3['u'],K4[’u']); 

V:=V+KUTTA(KI  [v ']  ,K2  [V]  ,K3[  Vj.K4  [  Vj); 

W:= W+KUTTA(K1  [’ w  ’  ]  ,K2  [Nv ' 1  .K3[\v ’]  ,Iv4  [’w’j ); 
X:=X-i-IvL’TTA(Kl[’x’],K2[,x’])K3[,x']JK4f’x']); 

Y:=Y-rIvL’TTA(Kl['y’J,Iv2[  y'J,K3[’y'],K-4|  y']); 

AU[I]:=U: 

AV[I]:=V; 

AW[I]:=W; 

AX[I|:=X; 

END;  {End  of  Outside  Loop} 

PRINT(AX, 250,1, 'S3’, C[3|); 

FOR  I:=l  to  N  DO 
BEGIN 

FOR  J:=l  to  100  DO 
BEGIN 

RUNGE(U, V,\V,X,Y, KAPPA, Kl); 

RUNGE(U+0.5*Kl[’u’],V-r0.5*Kl[’v’],W+0.5*Kl['w’],X+0.5*Kl['x'j. 

Y+0.5*Kl[V].KAPPA,K2); 

RUNGE(Up0.5*K2[’u’],V+0.5*K2[V],W+0.5*K2[\v’].X-i-0.5*K2['x’\ 

Y +0.5*K2[  V]  , KAPPA, K3); 

RUNGE  (U +K3  ( ’  u '  ] ,  V+K3  [  V  ] ,  W+K3  [  ’  w  ’  ]  .X-1-K3  [’  x  ’  ] ,  Y +K3  [  V  ] ,  IXAP  P  A .  K4 ) : 
U:=U+KL:TTA(Kl[V],K2(’u’],K3[V],K4[’u’]); 

V : = V  -fKU  TT  A(Kl  [  V]  ,K2  [  V  ]  ,I\3  [  V  ]  ,K4  [  V  ] ) ; 

W:= W+KUTTA(Kl  ( V]  ,K2  [  V]  ,K3  [\v’J  ,K4  [  V]) ; 
X:=X+KX'TTA(Kl[,x’],K2[,x’],K3[,x’],K4[’x’]); 

Y:= Y +KUTTA(K1  [V]  ,K2  [  V]  ,K3[  V]  ,K4  [  V] ); 

END;  {End  of  Inside  Loop} 

END;  {End  of  Outside  Loop} 

FOR  I:=l  to  250  DO 
BEGIN 

RUNGE(U.V.W.X.Y.IX-VPPA.Kl): 

RUNGE(UP0.5*Kl[’u'],V-r0.5*Kl[V},\V-f0.5*Kl[\v';  X-f-0.5*Kf'x';. 

Y+0.5*I\1  [  V]. KAPPA. K2): 

RUNGE(U+0.5:rK2[’u'j,V-f0.5*Iv2[’v’J.\V-f0.5*K2|'w'].X-f  0.5*K2,-x'jL 
Y+0.5*K2['y'],KAPPA,K3); 

RUNGE(L’+K3['u’],YpK3[’v'],W4-K3['\v'],X+K3[‘x']  A’-fK3[’y' j. KAPPA. Kl): 
U:=U-fKUTTA(Kl  [’u!J.Iv2['u-].K3['u’J.Kl[’u'j); 

Y:=Y+KUTTA(Kl  [V].K2['v'].K3[V].K4[V]): 


w : = W-rKUTTA(Kl  [\v'].K2[\v'].K3[\v'l  .K4[Nv'l ); 
X: =XtIvUT  T A(Kl  [  ‘  x  ’  ]  ,K2  [ '  x ;  j  ,K3  [ '  x '  ] .  K-4  [ '  x '  ’ ) :  ' 
Y:= Y-i-KL‘TTA(Kl  Py ']  .K2  [  V  ]  ,K3  f'y '  j  ,K4  ['y  ’] ) ; 
AU[I]:=U: 

AV[lj:=V; 

AW[Ij:=W; 

AX[I]:=X; 

END;  {End  of  Outside  Loop} 

PRINT(AX,250,1,’S4’,C[4]); 

PRINT(C,4T2,'LnA’.M[LOLA]); 

END; 

PRINT(M, LOLA, S, ’Alpha  vs.  omega', C [5]); 

TEXTMODE  (2); 

END. 


II.  Calculation  of  Steady  State  Values  from  Section  IV.B. 

Program  STEDSTAT: 

{$U+} 

{$1  GRAPH.P} 

LABEL  1,2; 

CONST 

al0=1.25E-f-06; 

a21=lE+04; 

a32=lE+06; 

omegal=2E-i-07; 

omega2=1000; 

omega3=lE-r07; 

dl=-2E+07; 

d2=2E+07; 

TYPE 

complex=array  [1..2]  of  real; 
buff=array[0..-406j  OF  BYTE; 
number=array(l..250]  of  real; 
gob=array[l..250]  of  integer; 
title=string[20] ; 

VAR 

BUFFER:buff; 

SAMPLE:byte; 

Wl,W2,W3,PMIN,PPLUS,\V,P0,Pl,P2.P3:number; 

T  rtitle; 

J,L.rlag:integer; 

Nl,N2,N3,N4,DET, 

ZEll,ZEr2,ZEl3,ZE21,ZE22,ZE23,ZE31,ZE32,ZE33:complex: 
M,DD,  DELTA, 

Yl  1  ,Y1 2,Y1 3  ,Y2 1  ,Y22  ,Y23,Y3 1  ,Y32,Y33 . 

AL 1 1 .  AL 1 2  ,AL  1 3  ,AL2 1  ,AL22  AL23  ,AL3 1 .  AL32 .  AL33 . 

D:real; 

I.ITl,IT2,IT3.IT4.IT5.IT6.IT7.ITS:complex: 

RTl.RT2.RT3.RT  l.RT5RT6.RT7RTS:real: 

PROCEDURE  r  key  in  (x:REALA'AR  z:COMPLEX); 

BEGIN 


i 


z[2]:=0.0; 

END; 


I 


I 


§ 


PROCEDURE  ikeyin  (x:REAL;VAR  zrCOMPLEX); 

BEGIN 

z[l]:=0.0; 

z[2]:=x; 

END; 

PROCEDURE  keyin  (zl:COMPLEX;  VAR  z2:COMPLEX); 

BEGIN 
z2[l]:=zl[l]; 
z2  [2]  :=  zl  [2] ; 

END; 


PROCEDURE  negate(y:COMPLEX;VAR  z:COMPLEX); 
BEGIN 

z[i];=-y[i]; 

z[2]:=-y[2]; 

END: 


PROCEDURE  con j ug at e(y : C OMP LEX; VAR  z:COMPLEX); 
BEGIN 

z[i]:=y[i]; 

z[2]:=-y[2]; 

END; 


P  FUNCTION  modulus(y:COMPLEX):REAL; 

[  BEGIN 

I  modulus;=SQR(y[l])-i-SOR(y[2]); 

J  A  END; 

I 

j  PROCEDURE  i n v e r t ( y : C O MP LEX ; V AR  z:COMPLEX); 

VAR  mag:REAL: 

BEGIN 

mag:=moduhis(y): 
z[l]:=y[l]/mag; 
z[2j  :=-y  [2]/ mag; 

END: 

PROCEDURE  add(yl.y2:COMPLEX:VAR  z:COMPLEX): 
BEGIN 

z[l]:= yi[i]+y2[i]; 


ttw  *ScS  ® 


z[2j:=yl[2]+y2[2}; 

END; 


I 


I 


e 


PROCEDURE  subtract(yl,y2:CO\IPLEX;VAR  zrCOMPLEX); 
BEGIN 

z[l]:=yl[l]-y2[l]; 

2  [2]  :=y  1  [2]-y  2  [2j ; 

END; 

PROCEDURE  m ul t ip ly (y  1  ,y 2 : C OMPLEX;VAR  z:COMPLEX); 
BEGIN 

z[l]:=yl[l]*y2[l]-yl[2]*y2[2]; 
z[2]  :=y  1  [1  j  *y2  [2]  +y  1  [2]*y2  [l] ; 

END; 

PROCEDURE  divide(yl,y2:COMPLEX;VAR  z:COMPLEX); 
VAR  mag:REAL;w  1  ,w2:C OMPLEX; 

BEGIN 

mag:=modulus(y2); 

conjugate(y2.wl); 

multiply(yl,wl,w2); 

z[l]:=w2[l]/mag; 

z[2]:=w2[2]/mag; 

END; 

PROCEDL-RE  scalarmult(x:real:y:complex;YAR  z:CO\lPLEX): 
BEGIN 
z[l]:=x*y[l]; 
z[2j:=x*y[2]; 

END; 

PROCEDURE  MINIMUXI(num:number;  VAR  min:real); 
BEGIN 
min:=num[l]; 

FOR  J:=l  to  250  DO 
BEGIN 

IF  num[J]  <  min  THEN  min:=num[J]: 

END; 

END; 

PROCEDURE  \L\XI\[UM(num:number;  VAR  max:real): 
BEGIN 

max:=num[lj; 

FOR  J :  =  1  to  250  DO 
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BEGIN 

IF  num[J]  >  max  THEN  max:=num[J]: 

END: 

END: 

PROCEDURE  SCALE(min.max,dx:real;  nummumber:  VAR  mark:gob); 
VAR  s:real; 

BEGIN 

s:=(max-min)/l30: 

FOR  J:=l  to  250  DO 
BEGIN 

mark[J]:=round((num[J]-min)/s); 

END; 

END; 

PROCEDLTIE  PRINT(num:number;  dx:real;  T:title;  VAR  par:real); 
VAR  min,max:real;  mark:gob; 

BEGIN 

GRAPHMODE; 

DRAW(  10, 190,260, 190,3); 

DRA\V(10,190,10,10,3); 

DRAW(260, 190.260, 10, 3); 

DRAW(10, 10,260, 10,3); 

FOR  J:=l  TO  9  DO 
BEGIN 
L:=190-J*13; 

DRAW(10.L,13,L.3); 

END; 

FOR  J:=l  TO  9  DO 
BEGIN 
L:=10+J*2o; 

DRA\V(L,190,L.1S7,3); 

END; 


MINIMLAI(num.min): 

\LAX3ML'M(num.max): 

SC.ALE(mi  n.  max.  dx.num.  mark): 
GOTOXY(l,25); 

\VRITE(min:10,'  !,max:10,’  ’,(max-min):10): 

GOTOXY(l,l): 

\\rRITE(((m  ax-*-min)/2):10,’  '.T): 
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FOR  J:=l  TO  249  DO 
BEGIN 

D  R.\\V(10tJ.  190-mark  rJ],ll-J.190-mark[J—l'. 3): 
END: 


READ(M); 

IF  M=1  THEN 
BEGIN 

FOR  J:=l  to  SO  DO 
BEGIN 

GETPIC(BUFFER.4*(J-l),0,4*(J-l).200); 

\VRITE(LST.#27/@’,#27,T'): 

WRITE(LST,#27,#75.#201,#0); 

FOR  L:=l  to  201  DO 
BEGIN 

SAMPLE:=BUFFER;(L-l)-r6j; 

WRITE(LST,CHR(S.\MPLE)); 

END; 

END: 


WRITELN(LST,#11.#11.#11.#11,#11.#10.#10, 

END; 

END; 


#10.#10.#27.'(a-): 


BEGIN 

DD:=(d2-dl)/250; 

DELTA:=dl: 


WRITELN(LST/dl=\dl/  d2=’.d2,’  deldel=\DD); 

WRITELN(LST,'Omegal=',omegal,’  Omega2=’.omega2,’  Omega3=’.omoga3  : 
\VRITELN(LST.’A10=‘.al0,'  A21  =  '.a21,'  A32=\a32): 

\VRITELN(LST,#11): 


IKE^ 

'IN(l.I); 

% 

RT1 

—omega 

il  *omegal : 

RT2 

=omega 

2*omega2: 

RT3 

=omeg? 

3*omega3: 

(Q 

N 1  [  1 

:=a2l; 

C 

N2[l 

:  =  a32: 

N3[l 

:=al0: 

Nil 

:  = :  1 10- 

a32; 

5 
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FOR  J:=l  to  250  DO 
BEGIN 


DELTA:=DELTA-DD: 


Nl[2]:=2*DELTA: 

N2[2]  :=2*DELTA: 

N3[2]:=2*DELTA; 

X4[2]:=2*DELTA: 

SCALAR\XXT(RT2*RT2.N2.ITl); 
MULTIPLY(IT  1  ,N3.ITS): 

{first  term  in  DET) 
X IUL  T  IP  L  Y(  Xl  ,X2 , 1 T 1 ) ; 
XXXTIPLY(ITl,N3JT2); 

SC  ALARXXXT  (RT 1  ,X2,IT  1 ): 
ADD(IT1,IT2,IT3); 

SCALARXXXT(RT3,X3,IT1); 

ADD(IT1,IT3,IT2); 

SCALARMULT(RT2*al0,IT2.IT7); 


ADD(IT8,IT7,IT6); 

{1st  and  2nd  terms  in  DET} 
NfLXTIPLY(X2,X3.ITl); 

MXXTIPLY(ITl,X4.IT2); 

S  CALARX  XX  T  (RT  1  .X3.IT1 ); 

.ADD(IT1  JT2.IT3); 

S  CALARX  XX  T  (RT3,  X2.IT  1 ); 

ADD(IT1.IT3,IT2); 


S  CALARX  XX  T  (RT2  *  a32.IT2.IT7 ): 


ADD(IT6.IT7.ITS): 

{1st,  2nd.  and  3rd  terms  in  E 
RKEYIX(RT  1-RT3.IT  1): 

X  XXT  IP  I  Y'(  X2 .  X-l ,  I T2) ; 

ADD(IT1,IT2,IT3); 

S  C  ALAR  X  XXT  (R  T 1  *  a  1 0  *  a32 ,1 T3  .IT7 ): 


ADD(IT3.IT7.IT6): 

{1st,  2nd.  3rd.  and  4th  terms  in  DET} 
RKEYIN(RT3-RTl.ITl): 

MULTIPLY(X3,M.IT2): 

ADD(ITl.IT2.IT3); 

SCALARMULT(RT3*alO*a3-2.IT3.IT7); 

ADD(IT6,IT7.ITS): 

{1st, 2nd. 3rd, 4th.  and  5th  terms  in  DET) 
\  lUL  T  IP  L  Y  ( X  2 ,  X  3 , 1 T 1 ) ; 

MlXTIPLY(ITl.N4.IT2): 

SCALARNILXT  (RT 1  ,X3,IT1 ); 

ADD(IT1,IT2.IT3); 

SCALAR.\1ULT(RT3,X2.IT2); 

ADD(IT3,IT2,IT1): 

SCALAR\fLXT(alO*a32,Xl,IT2): 

MEXTIPLY(ITl,IT2.IT3): 

ADD(ITS.IT3,DET);  {DET} 

RKEYIX(RTl-RT3,ITl); 

NIULTIPLY(N2,N4,IT2); 

ADD(IT1,IT2.IT3); 

SC\LARMLXT(RTl*a32,IT3sITS); 

RKEYIX(RT3-L_,ITl); 

MLXTIPLY(X3,X4,IT2); 

ADD(IT1,IT?,IT3); 

SCALARNrtXT(RT3*a32,IT3,IT7); 

ADD(ITS,IT7.IT6); 

MLXTIPLY(X1.X2.IT1); 

MlXTIPLY(ITl.X3.IT2); 

S  CALARX  H  XT  (R  Tl.X2.IT3); 

.ADD(IT2.IT3.rn ): 

S  C’ALARMl XT  (RT3.X3.IT2); 

,ADD(IT  1.IT2.IT3); 

S  C  ALAR  N  fULT  (RT2  .IT3.IT7): 


ADD(IT6.IT7.ITS): 


MULTIPLY(X2.X3,ITl): 

MLXTIPLY(ITLX4.IT2): 

S  C  ALARY  IULT(RT  1  .X3  .IT  1 ); 
ADD(ITl  .IT2.IT3): 
3CALARMULT(RT3.X2,IT2); 
ADD(IT2,IT3,IT1); 
SCALARMLXT(a32,Xl,IT2); 
MULTIPLYfITl  .IT2.IT7); 

ADD(ITS.IT7JT6); 

SCALARMULT(omegalXITl); 
MULTIPLYfITl, IT6.IT2); 
DI\TDE(IT2.DET,ZEll); 

Yl  1  :=ZE  1 1  [2] ; 

RKEYIX(RTI-RT3,IT1); 

MULTIPLY(X2,X4,IT2); 

ADD(IT1?IT2,IT3); 

SCALARMULT(a32,IT3.ITl): 

SCALARMULT(RT2,X2,IT2); 

ADD(ITl,IT2fIT3); 

SCALAR.\IULT(omegal*RT2.I.ITl); 

MULTIPLY(IT3.ITl,IT2); 

DrVTDE(IT2.DET,ZEl2); 

Y12:=ZE12[2]; 


{ZEll} 


{ZE12} 


ADD(X2.X3.ITl); 

MULTIPLY(I.ITl.IT2); 
SCALAR\XXT(omegal*RT2*RT3.IT2,ITl); 
DIYIDEfITl  .DET.ZE13); 

Yl  3 : = ZE 1 3  [2] ;  *  {ZE13} 

S  G  ALARY  ILXT  (RT2 .  N2  .ITS); 

RKEYIX(RTl-RT3.ITl); 

MULTIPLY(X2.Xl.IT2): 

.ADD(ITl  .IT2.IT3): 

S  C  AL  A  R  MIX  T  ( a  3  2 , 1 T  3 . 1 T 1 ) : 
ADD(IT1.ITS,IT7); 

MULTIPLYfI.IT7.IT8): 

SCALAlL\RXT(RTl‘o:n?ga2.ITS.IT7): 


DI\TDE(IT7.DET.ZE2l): 

Y21:=ZE2l[2]; 

MULTIPLY(X2,X4.ITl); 

SC\I_\RMLXT(a32.ITl.IT2); 

RKE'\lX(RTl'*a32.ITl): 

ADD(ITl.IT2.IT3); 

SCAL-\RMULT  (RT2.X2.IT2); 
ADD(IT3.IT2,ITl); 

S  C ALARN  ILXT  (RT2.ITl.ITS); 

SCAUXPAtLXT(a32.X4.ITl); 

RKEYIX(RT2.IT2); 

ADD(IT1.IT2.IT3); 

SC.YLYRMULT(RT3*alO.IT3.IT7): 

ADD(ITS,IT7JTS): 

RKEYIX(RTl*a32.ITl); 

SCALaRNIULT(RT2.X2.IT2): 

ADD(IT1,IT2.IT3); 

MULTIPLY(X2,X4,IT2): 

SCALARMLXT(a32.IT2.ITl): 

ADD(ITl.IT3,IT2): 

SC-AJL\RMULT(alO,Xl.ITl): 

MULTIPLY(IT1.IT2.IT7): 


ADD(IT6.IT7.IT$): 


S  C  ALARN I  LX  T  (o  m  eg  a2 . 1 . 1 T 1 ) : 
MULTIPLYiITS.it  1.IT2): 
DrUDElIT2.DET.ZE22): 
Y22:=ZE22’2’: 

RKEYIX(RT3-RT1.IT1): 
MLXTIPLYiX  1  .X2.IT2): 
ADD(ITl.IT2.IT3): 

S  C'ALAR  MLXT  (a  1 0. 1 T3 . 1 T  2 ): 

S  CAL  ARM!'  LT  i  R  T  2 .  X2  .IT  1 ) : 
ADD(ITl  .IT2.IT3): 

S  C  .ALAR  MLX  T  ■  o  m  c-  g  a  2  *  R  T  3 . 1 . 1 T 
\ ILXT I P  L  Y(  I T3 . 1 T  2  .ITT): 
DD1DE(IT1  .DET.ZE23): 
Y23:=ZE23r2;: 


{ZE21} 


{ZE22 } 


{ZE23} 


*  :y  • 


ADD(X2.X3.ITl); 

SCALARMLXT(-RTl*RT2*omega3.I.IT2^; 

MULTIPLY(ITl.IT2.IT3); 

Dr\TDE(IT3,DET.ZE3l);/ 

Y31:=ZE3l[2];  {ZE3l} 

RKEYIX(RT3-RTLIT1); 

MULTEPLY(Nl,N2.IT2); 

ADD(IT1,IT2JT3); 

SCALARMULT(alO,IT3,IT2); 

S  C  ALARX  IULT  (RT2  ,X2  .IT  1 ) ; 

ADD(IT1,IT2.IT3); 

SCALARMULT(RT2*omega3.LIT2); 

MULTIPLY(IT3.IT2.IT1); 

Dr\TDE(ITLDET,ZE32); 

Y32:=ZE32[2];  {ZE32} 

RKEYIX(RT  1-RT3.IT  1 ); 
\IULTIPLY(X2,X4JT2); 

ADD(IT1,IT2,IT3); 

SCAIAR\fULT(RTl*alO.IT3,ITS); 

RKEYIX(RT3-RTl.ITl); 

■  \IULTIPLY(X3,X4.IT2); 

ADD(IT1,IT2.IT3): 

SCALARMULT(RT3*alO,IT3.IT7); 

ADD(IT8,IT7,IT6); 


\  fULTEPLY(  X2  ,X3  ,IT  1 ); 
MULTIPLY(IT1,X4,IT2); 
SCAIARXIULT(RT1,X3,IT3); 
ADD(IT2,IT3.IT1); 

S  CALARm  u  It  (RT3  ,X2  JT2); 
ADD(IT1  .IT2.IT3); 

S  C  AL  AR  X  IULT  (R  T  2 . 1 T  3 . 1 T7 ) ; 


'f 

'♦p 


SCA1AR.\1LXT(RT3..\Y.IT2): 

ADD(IT3.IT2.ITl): 

\[ULTIPLY(ITl.Xl.IT2): 

SCAJL\RMULT(alO.IT2.IT7): 

ADD(IT8,IT7.IT6); 

\fLXTIPLY(IJT6JT7): 

S  CA1ARN  1XX  T  (o  m  e  g  a3 . 1 T  7 . 1 T  S ) : 

Dr\TDE(ITS,DET.ZE33); 

Y33:=ZE33[2]; 


{ZE33} 


ALII 

All2 

AL13 


.4121 

AL22 

.4123 


.4131 

.4132 

.4133 


=-alO/2-fa21/4-2*omegal*Yll—  omega2*Y2l; 

=al0-4  a21/2-2*omegal  *Yl24-omega2*Y22: 
=alO/2-a2l/4-2*omegal*Yl3-roinega2*Y23; 

=(alO-2*a2l4-a32)/4-2*omega2*Y21— omegal*Yll-rcmega3  *Y31 
=(-al04-2*a2l4-a32)/2-2*omega2*Y22-romegal*Yl24-omega3  ‘Y: 
=(-al04-2*a2l4-3*a32)/4-2*oinega2*Y234-omegal*Yl3— omega3’ 

=-a32/24-a21/4-2*omega3*Y3l4-omega2*Y2l; 

=-a324-a21/2-2*omega3*Y324-omega2*Y22; 

=-3*a32/2-a2l/4-2*omega3*Y334-omega2*Y23; 


D : = AL 1 1  *  (A122*  A133-AL23*  A132) 

-4112  *  (.4123*  .4131-.412 1  *  A133) 

4- .41 1 3  *  (.412 1  *  .4132-.4122  *  AL3 1 ); 

\Vl[J]:=((alO-a2l/2)*(.4l22*.4l33-.4l23*.4132) 

4-.4ll2*(A123*(a32-a2l/2)-(a21-al0/2-a32/2)*A133) 

4-.4ll3*((a21-al0/2-a32/2)*A132-Al22*(a32-a21/2)))/(2'D): 

V2[J’:=(  ALl  1  *((a2l-al0/2-a32/2)*A133-A123*(a32-a2l/2)) 

-(a  1 0-  a  2 1  /2  )*  (.-4123  *  A13 1-.412 1  *.4133 ) 

+ .41 . 1 3  *  (AL2 1  *  (a32-a2 1  /2)-(a21-alO/2-a32/2)*AL3 1  ))/(2  *D); 

\\’3[J;:  =  (ALl  l*(AL22*(a32-a2l/2)-(a21-alO/2-a32/2)*A.L32) 

4- .41 1 2  '(( a  2 1  -  a  1 0/  2-  a32  /  2 )  *  AL3 1  -.412 1  *  ( a32*  a2 1  /  2 ) ) 

-(a  1 0-  a2 1  /2)  *  (.411 1  *  AL32- AL22  *  AL3 1  ))/(2  *D): 

PMIX[J-:  =  (l-\Vl  [Jj-2*\V2[J]-\V3[J])/2; 

PPLUS  [.I] :  —  ( 1  -fW  1  [.J]  4-2  *\V2  [  J  ]  4~\V3  [  J  j  )/2; 


\V[J]:=PPLUS[Jj-PMIN[j;: 

P0[J]:=(l-3*\Vl[j;-2*\V2:Jj-W3:Jj)/4: 

P 1  [  J  j :  =  ( 1  ~\V1  [  Jj-2*\V2 r  J|-\V3  [  J] )/  4 ; 

P  2  [  J] : = ( 1  -f\  V  1  [  J  j  — 2  *  \  V2  [  J|-\ V3  [  Jj )/  4 : 
P3[J]:=(l-r\Vl[J]-r2*W2[J]-i-3*\V3;Jj)/4: 


END; 

PRINT(W1,DD.'\V1  vs.  Delta’. M): 
PRINT (W2.DD/W2  vs.  Delta’ .\I); 
PRINT(\V3.DD.’\V3  vs.  Delta’.M): 

PRINT(P\1IN.DD.'P minus  vs.  Delta’, M): 
PRINT(PPLUS.DD.'Pplus  vs.  Delta’, M); 
PRINT(W,DD.’W  vs.  Delta' .M); 

PRINT(PO.DD.'PO  vs.  Delta’, M); 
PRINT(Pl,DD,'Pl  vs.  Delta’.M); 
PRINT(P2,DD.’P2  vs.  Delta’.M); 
PRINT(P3.DD,'P3  vs.  Delta’, M); 


TEXTM0DE(2); 
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HI.  Calculation  of  Adiabatic  Variables  from  Section  IV. B. 3. 

Program  AD  LAB: 

{SU-} 

{$1  GRAPH.P} 

LABEL  1.2: 


CONST 

al0=lE-06: 

a21=l: 

a32=1.25E— 06: 
omegal  =  lE— 07 : 
omega2=lE  —  03; 
omega3=2E-r07: 
dl=-2E-r07 : 
d2=2E— 07 ; 

TYPE 

complex=array  [1..2J  of  real; 
buff=array[0..406]  OF  BYTE; 
number=array[l..250]  of  real: 
gob=array[l..250j  of  integer; 
title=string[20j; 


ALAR 

BUFFER:buff; 

S.AMPLE:byte: 

GAN  IX  LA  1 .  G.AN  1 1  GAN  1 3 .  G.AA  IN  LA2 .  G.AM 2  GAN  13 .  CLAN  IN  LA3  ,E  T  A 1 .  E  T  A  2 .  E  T  A ; 
DELT.Al.DELTA2.DELTA3dvIl.XI2.XI3.PSI. PHI. RP.RM. 

PP,  PM,  TEST  :n  umber; 

T :  t  i  1 1  e : 


J.L.fiag:integer; 

XLX2.X3.Xf.DET. 


R 1 1  .R 1 3  .R  M  ,R2 1  .R22  .R23.R2-1  .R25  ,R3 1  .R32.R33.R3-1.R  3-5 . 


R-ll.R-12.R-13,Rd-l.R-15:complex: 


M.DD. DELTA. 

AL 1 0  ..AL 1 1 .  AL 1 2 .  AL 1 3 .  AL20 .  AL2 1 .  AI  .22  .A  L23 . 
B11.BI2.B21. B22.BXT0.CT1. C12.C13.C20.C21.C22.C23. 
D:real: 


UT1,IT2.IT3,  IT-1,  IT5,IT6.IT7.ITS:complex: 
RT 1. RT2. RTS. RT-1.RT5.RT6.RT7. RTS:  real: 


PROCEDURE  r key in  ^x:REAL:YAR  -COMPLEX); 
BEGIN 


z!l,:=x: 


z:2;:=0.0: 


ENT): 


PROCEDL'RE  ikeyin  (x:REAL:VAR  z:COMPLEX); 
BEGIN 
z[l];=0.0: 
z[2]:=x: 

END: 


PROCEDL'RE  kevin  (zl:COMPLEX:  VAR  z2:COMPLEX) 


BEGIN 

z2[l]:=zi;ij: 


z2f2]  := z  1  r2l ; 


l-J 

END; 


PROCEDURE  negate(y:COMPLEX;YAR  z:COMPLEX): 
BEGIN 

z[i]’-=-y[i]; 

z[2l:=-y:2l; 

END: 


PROCEDURE  conjugate(y:COMPLEX:YAR  z:C’OMPLEX); 
BEGIN 


z[ll:=v[l’; 

l  -  l  i 


END: 


FUNCTION  modulus(y:COMPLEX):REAL; 
BEGIN 

modulus : =S  OR  (y[lj)-f-SQR(y[2 1 ); 

END: 


PROCEDL'RE  invert(v:COMPLEX:\'AR  z:COMFLEX) 
WAR  mag :REAL: 

BEGIN 

mag:=moduIus(y); 
z[l]:=y[lj/ mag; 
z[2]:=-y[2j/mag: 

END: 


PROCEDURE  add(vl.y2:COMPLEX:YAR  z:CO\IPLEX  : 
BEGIN 

z[l’:=yE  1  ■— y2;l’: 
z;2!:=yl  2l—y2  2  : 

END: 

PROCEDURE  subtract(yl.y2:COMPLEX:VAR  z:COMPLEX,: 
BEGIN 

z  [  1  ]  :=y  1  [  1  ’ -y  2  ( 1  ]: 
z[2]:=yl  ;2l-y2[2j: 

END: 

PROCEDURE  multi ply(yl.y2:COMPLEX:\".XR  z:COMPLEX;: 
BEGIN 

z  [  1  j  :=y  1  [  1  ]  *y2  [  1  j-v  1  [2]  *  y2  [2] : 
z[2':=yl[l]*y2[2j-ryl[2]*y2[l]; 

END: 

PROCEDEDE  divide(yl,y2:CO\IPLEX;YAR  z:C'OMPLEX): 
VAR  mag:REAL:wl,w2:CO\IPLEX; 

BEGIN 

mag:=modulus(y2); 
conjugate(y2,\vl): 
mu!tiply(yl  .w  1  ,w2); 
z[lj:=w2[l]/mag: 
z[2j:=\v2[2]/mag: 

END: 

PROCEDURE  seal  annul  t(x:reaI:y:comp!ex:\'.-VR  z:CO\IPLEX): 
BEGIN 
z[l]:=x*y;l]: 
z[2]:=x*y[2;: 

END: 

PROCEDURE  M I NI M  LX  I  ( n  u  m :  n  u  rn  b  e  r :  VAR  mimreal): 
BEGIN 
min:=num[lj: 

FOR  J:  =  1  to  250  DO 
BEGIN 

IE  num[J]  <  min  THEN  mir.:  — runr/.I’: 

END: 

END; 


PROCEDURE  MAXIMUM(num:nuinU'n  \'AR.  max:real) 


BEGIN 

max:=num;l’; 

FOR  J:  =  l  to  250  DO 
BEGIN 

IF  num[J|  >  max  THEN  max:=num;J|; 

END; 

END; 

PROCEDURE  SCALE(min.rnax.dx:real:  num:number;  VAR  marlcgobj 
VAR  s:real: 

BEGIN 

s:=(max-min)/lS0: 

IF  s=0  THEN  s:=l; 

FOR  J:=l  to  250  DO 
BEGIN 

mark[J]:=round((num[J]-min)/s); 

END; 

END; 

PROCEDLRE  PRINT(num:number;  dx:real;  T : t i 1 1  e ;  VAR  par;rsal): 
MAR.  min,max:real;  mark:gob: 

BEGIN 

GRAPHMODE; 

DRA\V(10, 190.260.190.3); 

DR\\V(10, 190.10, 10,3); 

DR\W(260. 190.260, 10,3): 

DRAW(lO.  10,260, 10,3); 

FOR  J;=l  TO  9  DO 
BEGIN 
L:  =  190-J*1S; 

DRA\V(10.L.13.L,3); 

END; 

FOR  J:=l  TO  9  DO 
BEGIN 
L:=10+J*25: 

D  R  A  \  V  ( L .  1 9  0 . 1, .  1 S  7 . 3  j : 

END; 

MINTMUMfaum.inin): 

NLAXIMI.’Mfnmn.inax); 

SCALE!  in  i  n .  rn  ax .  d  x ,  n  u  in  .mark): 

GOTOXY(  1 .25): 


\VRITE(min:lO.’  \max:10.'  ’.(max-minLlO): 
GOTOXYi  1.1): 

\VRITEiT(  max— m:a)/2):10.'  -T): 


FOR  J:=l  TO  249  DO 
BEGIN 

DRAW(10-f  J.190-mark[J].ll4-J.190-mark[J-r  l].3); 
END; 


READ(M); 

IF  M=1  THEN 
BEGIN 

FOR  J:=l  to  SO  DO 
BEGIN 

GETPIC(BUFFER,4*(J-l),0,4*(J-l).200); 

MPITE(LST.#27/@’.#27.T): 

\VRITE(LST.#27,#75,#201.#0); 

FOR  L:=l  to  201  DO 
BEGIN 

SAMPLE ; = BUF  F  ER  [(L- 1 )  • 4-6  j ; 
WRITE(LST.CHR(S  AMPLE)); 

END; 

END; 


\\PJTELN(LST.#11.#11.#11.#11.#11.#10.#10.#10.#10 

END; 

END; 

BEGIN 

DD:=(d2-dl)/250; 

DELTA:=dl; 

WRITELN(LST AD  IAB' ): 

WRITELN(LST,'dl=\dl.'  d2=\d2/  doldeI  =  ’.DD): 
WRITELN(LST,’0megal  =‘.cmc-gal Omega2=D.onifi::i2 
\  VR I T  E  L  N  ( L  S  T . '  A 1 0 = a  1 0 . '  A21=\n21/  A32=\a32): 

\VRITELN(LST.#ll); 

IKEYIN(l.I); 


RT1 

RT2 

RT3 


=omegal  ‘omcgal ; 
=omega2*omega2: 
= o  in  e  g  a3  *  o  n  i  ?  ■  g  a  3 : 
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Xl[lj:=a21: 

X2[l’:=a32: 

X3fl  j:  =  alO— a21: 

X4fl  f:=alO— a32: 

FOR  J:=l  to  250  DO 
BEGIX 

DELTA:— DELTA-DD; 

Xl[2j:=2*DELTA: 

X2  [2  j : = 2  *  D  EL  T  A : 
X3[2j:=2*DELTA; 
X4[2]:=2*DELTA: 

RKEYIXfomegal/alO.ITl): 

MLXTIPLY(ITl.I.Rll); 

RKEYIX(-omega2/(2*alO).ITl): 

MLXTIPLY(IT1.LR13): 

MLXTIPLY(R13,I.R14): 

N ILX  T  IP  L  Y  ( X  2 .  X  3 . 1 T 1 ) : 
MLXTIPLY(IT1.X4.IT2): 

SCALARN ILXT  (RT 1  .X3.IT3): 
ADD(IT2.IT3.ITl); 

S  C  AL A  R  MIX  T  ( R  T  3 .  X  2 . 1 T  2 ) : 
ADD(IT1.IT2.IT3): 

S  CLAL AR  MXX  T  ( a  3  2  -  a  2 1 . 1 T  3 .  ITl ) ; 

N  ILXT  IP  L  Y(  X2 .  X3 . 1 T  2 ) : 

S  C  A  L  AR  M  t  X  T  ( R  T  2 . 1 T  2 . 1 T  3 1 : 
ADD(IT1.IT3.DET): 

MLXTirLYfX2.X4.ITO; 

RKEYIXi  RTl  ,IT2  ): 

ADDfITl  .IT2.IT3): 

S CAI >AR .\ HXT (a32-a2l.IT3.rn  : 

5  CAI . AR \  ILL T  ( R T2 .  X2 . 1 T 2  : 
ADDfITl.  IT2.IT3): 

S  C  A I.  A  R  ML '  L  T  ( -  o  m f  g  a  2  /  2 .  IT  3 . 1 T  2 1 : 
MLXTIPLY(LIT2.m): 

D m I) E ( ITl . D i ;  r .R 2 1 ) : 


lor, 


V* 


KM 


*V*V 


S  C  ALAR  \  X ' X ' T I  o  m  e  5  a  2  *  R T  3 .  X  2 . 1 T  1 ' : 
MLXTIPLYiI.ITl.IT2): 
DATDEiTT2.DET.IT3): 
ADDiIT3.R21.R22!: 


S  C  ALAR\  fLX  T  ( - 1  .R2 1  ,R25 1: 


MULTIPLY(X2.X4.ITl): 

RKEPTX  (R  T 1  -  RT  3 . 1 T  2 ) ; 

.ADD(IT1.IT2.IT3): 

SC.\LAR.MLXT(a32-a21.IT3.ITl): 

S  C  ALAR  MUL  T  (R  T  2 .  X  2 . 1 T  2 ) : 

XDD(ITl.IT2.IT3): 

SC\L.ARMLXT(omegal/2.IT3.ITl): 

MLXTIPLY(I.IT1.IT2): 

DATDE(IT2.DET.R23): 


MULTIPLY(I.R23.R2-l): 


S  C  XL . AR  MLX  T  (-  R  T  2  *  o  m  e  g  a  3  /  2 .  X  2 . 1 T 1 ) : 

MXXTIPLY(I.ITl.IT2); 

DATDE(IT1.DET.R31): 


;C.ALARMLXT(- 1  .R31.R35): 


MLXTIPLY(X2.X3.IT1): 

MLXTIPLY(IT1.X-LIT2): 

S  CALARN  XXT  ( R 1'  3 .  X2  .IT  1) 
.ADD(ITl.IT2.IT3}: 
SCALARMIXT  (RT  1  .X3.IT2; 
ADDiIT3.rr2.ITR: 
SCALARMULTiomegii3.IT  1 .1 
MI.T/nPLY(  I.IT2.IT  1  j: 

D I  \T  D  F.  ( TIT .  I)  E  T .  I T  2 1 : 

AD  D  ( I T  2 .  R  3 1 .  R  32 ; : 


21: 


ADD(X2.X2.rn  j: 
SCAI.ARMM.T  a  1 
MULTIPLY; I.rr2.m  >: 

I)IYiI)E:riT.I)ET.R33i: 


::i«  o 


rr  1  IT" 


MI  'I/I  ILLY;  I.R33.R::  l  : 


l‘i; 


•:*,*  v;v,' 


PJ\EYIX(-omegal‘oinega2*ornega3:’‘(  a21-a32)  2.IT1 ' 
NILXTIPLYiI.ITl.IT2): 

DINTDE(IT2.DET.R4l): 

S  CALARN  fLXT(-l.R41  .R45) : 

SC.ALARNrLXT(-omegal*omega2xcmega3.X3.ITl): 

MLXTIPLY(I.IT1.IT2); 

Dr\TDE(IT2.DET?ITl): 

ADD(IT1.R41.R42); 

MLXTIPLY(X3.X4.ITl): 

RKEYIX(RT3-RTl.IT2); 

.ADD(ITl.IT2.IT3): 

SCALARNItXT((a32+a2l)/2.IT3.ITl); 

SCALAR  MLXT(RT2/2.X3.IT3): 

ADD(ITl.IT3.IT2): 

SCALAR\lLXT(-omega3,IT2.ITl): 

MLXTIPLY(I.ITl.IT2); 

DINTDE(IT2.DET.R43): 

MLXTIPLY(I.R43.R44); 


— al0-f-2*omegal  *Rl  1  [2j-omega2*R2I  [2?; 
=a2l/2-omega2*R22[2j; 
=2*omega3*R3l[2]-omega2*R21/2i: 
=a32Aa21/2A2*omega3‘R32)2)-omega2‘R22’2 


ALIO: 

.4X20 

ALII 

AL21 

.4X12 

.4X22 

AL13 

AL23 


=-(alO-a21  /  2 )  /  2 : 

=-(a32-a21/2)/2: 

=  ( a  1 0*  a2 1  /  2 — 2  *  omega2  *  R  25  [2]  )/2 : 

— -u32/2-f  a2 1  /4-2 ‘omegaS  *R35/2; -omega2  'R25  >2; 


=cmega2‘R23[2 
=omcga2kR23)2 
— G!neaa2*R24  [2 


-2*o:::ogal  ‘R13/2  : 
-2 ‘omeii'.o * R33  2  : 


Al.23:— om-'g:i2 * R2  1/2  -2 *<.nKga3 ‘  1  ’34  2  : 
B:  =  B1 1  rB22-Bl2“B21: 


CIO 

__i  i  >•>■> * 

AI.10-R1 

O  * 

\I 

20  CR 

.  i 

C20 

—  i.Bll* 

AL20-P>2 

1  r 

Al 

lOi/B 

Cl  1 

_ !  I  )-)0  * 

- \  I  >-- 

A  L 1 1  -  B 1 

■  )  * 

\1 

21  ;./B 

C21 

1  • 

A  1.2 1-P>2 

1  * 

Al 

1 1  :/B 

C 1 2 

'  1 522  * 

AL12-B1 

•)  » 

A! 

22  X. 

=iBll*AL22-B21“ALl2'  B 
=(B22*  AL 1 3-B 1 2* AL23 !,  B 
B 1 1  r  AL23-B2 1  \AL13 /  B 


GAG  IG  LAI  LJ’  :=a21  — 2*om?gal  ‘!R23’2;— R2L2;  ‘  C 12-R22  2'  ‘  C22 
-(omega2/alO')*(2*omegal*Cl2-omega2) 

-2*omega3*(R43[2i  — R4 1  2  * C 1 2— R42 , 2.  * C 2 2  i : 

G.4MM\2[.J;  :=a2 1-2‘omegal  *  (R24  [lj  -R2 1 : 1  ] « C 1 3-R22  Y*  C23  i 
J-(RT2/alO)-r2*cmega3*iR44ll R41  tl j * C 1 3 — R4 2 _  1  / C 23  : 

GAMNL\3[J::=a21*(l-2*C21;— 2*cmega2*(R25/2]— R2l’2^  *Cll— R22/2"  *  C21 

GAGUGAGI3[.Jj :  =  GAGENLA1 :  J;/GAG1G  IA3;  J ; : 

GAM2GAG  13  [  J  ]  := GAG  L\  LA2 :  J]  /  GAG  IG  LAS :  Jj : 


ETAl[Jj:=-omegal*(R21  [2]*C10-i-R22.2]*C20)— (omega  1  *omega2/al0  "CIO 
+omega3*(R4l[2]*Cl0—R42.2]‘C20); 

ETA2;J]:=omegal  *(R21  [lRClO— R22!l:*C20)-cmega3*(R41  Y  *ClO— R42  Y 

ETA3[J]:=a21*(C20-0.5)-omega2*(R2l[2]*Cl0-rR22;2;,C20): 

DELTAl[J;:=DELTA-omegal*(R24[2]-rR21 :2;*C13— R22’2/C23 1 

■f(omegal *omega2 /al 0) * C 1 3-fomega3 * (R44 [2]  — R4 1  [2< *03— R 4 2  2  *Ci 

DELTA2[  Jj  :=DELTA-omegal  *  (R23  [l ;  -r  R2 1 Y  *  C 1 2  4-R22  Y  *  C22  i 
-romegaS*  (R43ll '  — Rl  1 Y  *  C 1 2— R42 f  1 1  *  C22): 


DELT A3 : .] i : = a2 1  kC22-omega2'(R23f2' '.4-R21  '2'.  *C12-R22'2'  *C22  : 

Xn[Jj:=-omegal  *(R25;2>R21  Y'0 1 +R22X  *C21 ) 

Vfomegal  *om*ga2/alO)*C  1 1  AoniegaS  *(R  15[2; -i-Rl  1  Y  *  C 1 1  R  12Y '  C2 


XI2LJ :  Iomega  1  * (R25 [  1  ]  AR 2 1  [  1  j  *  C 1 1 -R22  Y  * C2 1 ') 
-omega3*(R  loY-R 4 1  Y  *C'  1 1  -*-R  12Y*C21 ): 


X13.J  a2 1  *C23-o:noga24 1R2  !  2]  —12 2 1  2’  *  C‘ ]  R22  2’  *  C23 

;o:=- ■1/(GA.\1GIA1;J>GAMMA2;J>4 'DELTA i;.i;*I)EI.-TA2;.r.it 


ii;  tagge 

4-RTo  *  (’delta  3 Ll ;  *  (0.5 '  G  AMM  A2  '.r * 


I  .T  \ l  M  -  D1  !  r A  !  A  \2 


vvv  NX* 


4IPIJII 


MRPKIPVFIMmnilW 


£ 


— XI3\J’*i  0.5  “GAMA  L\1  J' *ETA2’J’-DELTA2' 


■  c 


i  a: 


■a. 

A 


P5i;j;:=0..5-G.V.\L\L\3;j; 

-RT5*(  DELTAS' T  * .  0.5  *  CAMNLV2  T  *XIl  T  -DELTA1 ;  J‘  ‘XT2 

-xi3;j;*to.o*  gammagj; -xi2;j;-delta2;.]; -Xu 


RT4:=G.\AL\  LAI  ;.j;*GAM\L-V::.J’'G.\-\L\L\3;j;/S-Xl2;.r‘DELTAl'j' ‘DELTAS'.]' 
-XI 1  '-E  *  DELTAS  3J ■  *XI3 )  J’- GAA  IN  LAS :  Ji  ‘XI 1  ).L  ‘ DELTAS  ’  J' ,  2  ’ 

— D E L T A 1  J •  ‘ D E L T A 2  J ’  ‘ G .AN EN L\3 ■  J •  / 2 - G .AN IN L\ l’ .E x XI 2  •  J  ‘ XI 3  J ’  2: 


RM;j;:=o.5«iPsi;j;-PHrj;): 

RP;j;:=o.5‘!psi'j;-pin;j;;; 

pp;j’:=rp  ;.j;  ;  <.  rp;j;-rm;j;  • ; 
p\i;j;:=rm;j;/iRp  ;j;~rm  ;j  > 

test;j;:=pp;j;-pm;j;: 

END: 


P RIN’T ( GAN IN  LA  1  .DD .  ’  G.  VMN LA 1  VS.  Deita’.M/: 
PRINT! G AX IN  LAS .  D  D . '  G  AMN  LA 2  \'S.  Deita’.M): 
PPJNT(GAM.\LA3.DD.'GAM\LA3  VS.  Delta' .M): 

P R1NT( G.AN  1 1  G.AM3  .DD  .'GAA 1 1  / GAM 3  VS.  Delta’ .M): 
P  R I N  T  ( G  AN  1 2  G  .-AN  1 3 .  D  D . '  G  AM  2  /  G .AN  1 3  VS.  Deita'.M): 
P  RINT(DE  LT  A 1  .D  D .  ’D  -:•!  t  a  1  VS .’  Del  t  a’  .M): 

PRINTi  DELTAS. DD.'DeltaS  VS.  Deita’.M): 

PRINT!. PHI.DD.’PHI  VS.  Deita’.M): 

PRINT1. PSI.DD.'PSI  VS.  Deita'.M); 

PRINTi  RM.DD.'RM  \'S.  Deita'.M!: 

PRINT!, RP.DD.'RP  VS.  Deita'.M); 

PRINT(PM.DD.’PM  VS.  Deita’.M}: 

PRINTi  PP.DD.'PP  \‘S.  Deita'.M): 
PRINTiETAS.DD.'ETAS  \'S.  Deita'.M): 
PPTNTiDELTAS.DD. 'DELTAS  VS.  Deita'.M): 

PRINTi  TEST. DD. 'TEST  VS.  Deita'.M): 


TEXT  MODE  2): 
END. 
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